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I, Shankar Subramaniam, PhD., declare as follows: 

1) I am currently a Professor at the Departments of Bioengineering, Chemistry and 
Biochemistry at the University of California, San Diego. I was a Professor in the Departments 
of Physiology, Biochemistry, Chemical Engineering and Electrical and Computer Engineering at 
the University of Illinois prior to moving to UC San Diego. Previously, I have held Assistant and 
Associate Professorships at the University of Illinois. 

2) I obtained a Bachelors of Science majoring in Chemistry, Physics and 
Mathematics in 1972 from Osmania University, Master of Science in Chemistry in 1974 from 
Osmania University and a Ph.D. in Chemistry from the Indian Institute of Technology Kanpur in 
1 982. I have authored or co-authored numerous papers in the areas of genomics and 
bioinformatics. My curriculum vitae and list of publications is attached hereto as Exhibit 1 . 

3) I have reviewed sections pertinent to my below attestation of the above-identified 
patent application. Specifically, I understand that the application describes and claims computer 
methods for producing a genome specific stoichiometric matrix and computer methods for 
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producing an in silico representation of a microbe. The methods rely, in part, on the assignment 
of function to a gene or open reading frame based on nucleotide or amino acid sequence 
homology. 

4) I have read the Office Action mailed June 18, 2004, issued in connection with the 
above-identified application. I understand that the claims have been rejected, in part, because it 
is alleged that the application does not describe percent homology criteria for functional 
categorization of open reading frames based on sequence homology and that the application does 
not describe criteria for inclusion or exclusion of a gene or open reading frame as a metabolic 
gene based on sequence homology. While all open reading frames in a genome cannot be 
automatically assigned function, it was well established at the time the priority application was 
filed on February 2, 1999, that those open reading frames that have appropriate sequence 
homology with genes from other microbes whose function was previously assigned could be 
characterized to possess the same function. While the degree of sequence homology is arguably 
subjective, databases such as COG (Complete Groups of Orthologous Genes, National Institutes 
of Health Website that uses bi-directional BLAST; http://www.ncbi.nlm.nih.gov/COG/; Science 
1997 Oct 24;278(5338):631-37) have criteria that are routinely deployed for function 
identification through homology. 

5) My understanding from reading the application is that, at the time of the priority 
application, functional assignment through homology was a routine procedure for open reading 
frames that display homology, and it is a superfluous need for the procedure described in the 
application to be set forth in more detail than that currently specified. Further, once function 
identification is made it is common place to characterize the function as pertaining to metabolism 
from legacy knowledge of cellular metabolism. 

6) For example, the application describes that following identification of an open 
reading frame, well established algorithms can be used to determine the extent of similarity 
between a given sequence and a gene/protein sequence deposited in the worldwide genetic 
databases. These algorithms are described at page 7, second paragraph, and include the BLAST 
and FAST family of programs. 
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7) The application also describes at, for example, page 7, second paragraph, that if a 
coding region from a gene is homologous to a gene within one of the sequence databases, the 
open reading frame is assigned a function similar to the homologously matched gene, 

8) The application further describes at, for example, pages 13-14 in Example 1 ? that 
the genetic sequence and open reading frame identifications and assignments are readily 
available for a majority of the E. coli genome from a number of on-line locations, citing to both 
the TIGR and the E. coli Genome Project web sites for obtaining this information. Genome 
annotation data obtained from the latter site was used to obtain those genes involved in cellular 
metabolism. 

9) Based on the above descriptions, it is clear that for this aspect of the claims, all 
that is necessary with respect to assigning function to encoded proteins is to perform a sequence 
homology search against a gene or protein sequence database using well known algorithms or 
programs such as BLAST or FAST to identify if a coding region is homologous to a known gene 
within a sequence database. Identification of a homologous sequence assigns the protein 
function of the known homologous sequence to the protein encoded by the query coding region 
sequence. Accordingly, a positive sequence homology hit to a known sequence of known 
function results in the assignment of the known function to the query coding region. 

10) Based on the above descriptions, it also is clear that in order to practice this aspect 
of the invention as claimed, all that is necessary with respect to determining if the assigned 
function of an encoded protein is involved in cellular metabolism is to select those open reading 
frames that have been assigned a function involved in metabolism. As indicated at pages 13-14, 
the genome annotation data of the known gene will, in many cases, already contain such an 
assignment. For example, the annotation data can indicate that a particular gene is involved in 
glycolysis, one of many well known metabolic processes. In cases where the annotation of a 
known gene specifies only the name of the protein or a biochemical reaction, published literature 
can provide information on functional association of the gene with metabolism. 
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11) Further, at the time the priority application was filed, any experimentation 
necessary to obtain a sequence homology search and assign a function based on homology to a 
known gene was predictable and routine. For example, the BLAST and FAST family of 
sequence alignment algorithms were well known and accepted standards in the art. The percent 
homology results produced from these programs are based on alignment criteria selected by the 
user. It is also common to use an E-value that seeks if the two sequences are similar by random 
chance. It is routinely accepted in the protein research area, that sequences that are likely to be 
homologous by random chance by one in a tenth of a million are likely to be orthologous or 
paralogous sequences. Accordingly, the predictability of selection and assignment of function 
based on homology to a known gene is set by the user and can be based on well established 
criteria. 

I hereby declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true. 



Date: TXlC ZJJ ^ 0 ij 



Shankax Subramaniam, Ph.D. 



CURRICULUM VITAE (SHORT VERSION) 



SHANKAR SUBRAMANIAM 



RESEARCH PROFILE: 

Shankar Subramaniam is a Professor of Bioengineering, Chemistry and Biochemistry and 
Biology and Director of the Bioinformatics Graduate Program at the University of California at 
San Diego. He also has adjunct Professorships at the Salk Institute for Biological Studies and 
the San Diego Supercomputer Center. He is also a Guest Professor at the Center for Molecular 
Biology and Neuroscience at the University of Oslo in Norway and Professor at the Center for 
Cardiovascular Bioinformatics and Modeling at Johns Hopkins University. He is a fellow of the 
American Institute for Medical and Biological Engineering (AIMBE) and is a recipient of 
Smithsonian Foundation and Association of Laboratory Automation Awards and his research 
work is described below. Subramaniam was honored in 2002 as "High Performance 
Computing's Highest Guru" by Genome Technology Magazine. He is on the Editorial Board of 
the following Journals: Bioinformatics, Journal of Molecular Recognition and the journal, 
Molecular and Cellular Proteomics. He is on several NIH advisory committees and panels. His 
research spans diverse areas of bioinformatics and systems biology. 

RESEARCH INTERESTS OVERVIEW 

The sequencing of several mammalian genomes provides a basis for understanding the 
systemic functioning of living cells. The "omics" technologies have begun to produce vast 
amounts of context-specific biological data 

Thus, future developments in genomics, and the applications that derive from genomics, will be 
dependent upon the scientific progress at the interface of three major disciplines; biology, 
engineering, and computer science. My laboratory works in this interdisciplinary areas of 
Bioinformatics and Systems Biology. 

Bioinformatics characterizes the flow of information in living systems. This is schematically 
represented below. 



Specific projects our laboratory is associated with in this area include, 
Genome Annotation 

1. Design of a novel algorithm for generating protein family specific scoring matrices and 
applications of this scoring scheme in identifying orthologous and paralogous proteins in 
genomes. 

2. Design of nucleotide and protein profiles that will serve to develop a high throughput 
method to identifying genes in genome sequences. This has implications for expression 
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profiling technology using oligonucleotides for assessing cell function. 
3. Design and development of high throughput pipeline for annotating genomes. 

Protein Sequence-Structure Mapping 

1. Development of strategies for obtaining unique protein fragment structures that will serve 
as the basis set for modeling three-dimensional structures of all soluble proteins. 

2. Characterizing protein folds using fragment-based methods. 

3. Obtaining statistical- potentials (pairwise atomic PDFs) for defining energy landscape of 
proteins. 

4. Using global optimization strategies for modeling structures of proteins from their primary 
sequences. 

Functional Genomics 

1 . Development of new methods for gene expression analysis. 

2. Reconstruction and biochemical analysis of transcriptional networks. 

3. Functional analysis of transcription factors. 

Reconstruction and Modeling of Biochemical Pathways 

1 . Reconstruction of disease networks pertaining to insulin resistance. 

2. Design and development of methods for analyzing cell signaling pathways from 
expression profile and proteomic data. 

3. Identification of protein partners in cell signaling pathways. 

4. Computational modeling of pathways. 

Infrastructure for Biological Databases, Analysis Tools and Interfaces 

1. Extensions of the Biology Workbench (http://workbench.sdsc.edu ) developed in my 
laboratory. 

2. Design of novel tools for expression profile and proteomic analysis in the Workbench. 

3. Creation of automated pipelines for genome, gene and protein analysis. 

4. Design and development of Laboratory Information Management Systems. 

5. Development of the Information Management System for Alliance for Cellular Signaling and 
LipidMaps projects. 

http://qenome.ucsd.edu/ 

http://www.siqnalinq-qatewav.org/ 

http://www.lipidmaps.org/ 

http:// www.mitorproteome.org/ 

http://www.bioinformatics.ucsd.edu/ 



Funding Acknowledgement: NSF, NIH, State of California and Corporate 
AFFILIATIONS: 
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Professor 

Department of Bioengineering 
Department of Chemistry & Biochemistry 
San Diego Supercomputing Center 
Salk Institute 
Director 

Bioinformatics Program, UCSD 
Guest Professor 

Center for Molecular Biology and Neuroscience 
University of Oslo, Norway 

EDUCATION: 

B.Sc, (Chemistry) Osmania University 1972 (Recipient of gold medal for ranking first in the 
University) 

M.Sc, (Physical Organic Chemistry) Osmania University 1974 (First Class with distinction) 
Ph.D., (Chemistry) Indian Institute of Technology Kanpur 1982 



PROFESSIONAL EXPERIENCE: 



1999-present Professor, Departments of Bioengineering, Chemistry and Biochemistry, 

Director, Bioinformatics Graduate Program University of California San Diego 
and Senior Fellow, San Diego Super Computer Center. 

1998-1999 Professor, Departments of Biochemistry, Molecular and Integrative Physiology, 
Chemical, Electrical and Computer Engineering, Center for Biophysics and 
Computational Biology, Beckman Institute and National Center for 
Supercomputing Applications. University of Illinois at Urbana-Champaign. 

1996-1997 Associate Professor, Department of Molecular and Integrative Physiology, 

Department of Chemical Engineering, Center for Biophysics and Computational 
Biology, Beckman Institute and National Center for Supercomputing 
Applications. University of Illinois at Urbana-Champaign. 

1991-1996 Assistant Professor, Department of Molecular and Integrative Physiology, Center 
for Biophysics and Computational Biology, Beckman Institute and National 
Center for Supercomputing Applications. University of Illinois at Urbana- 
Champaign. 

1990-1990 Visiting Scientist, Princeton University and Senior Research Scientist, Squibb 

Institute for Medical Reseach. 
1986- 1989 Assistant Director for Scientific Development, IMD and Visiting Assistant 

Professor of Chemistry, University of Houston. Work carried out in the group of 

Professor J. Andrew McCammon 
1985-1986 Postdoctoral Research Associate, University of North Carolina, Chapel Hill. 

Work carried out in collaboration with Profs. Jan Hermans and Max Berkowitz. 
1984-1985 Postdoctoral Research Associate, University of North Carolina, Chapel Hill. 

Work carried out in the group of Prof. Robert G. Parr. 
1982-1984 Lecturer in Chemistry, St. Stephen's College, Delhi University. 
1979-1982 Senior Research Assistant, Indian Institute of Technology Kanpur. 
1976-1979 Research Assistant, Indian Institute of Technology Kanpur. 
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HONORS: 

Genome Technology All Star Award, 2002 

Association of LabAutomation Award, 2001 

Elected Fellow, Institute for Biomedical Engineering, 2000 

The College of Engineering, University of Illinois Advisor's List, 1997. 

Smithsonian Foundation Citation for Innovation in Computing, 1997. 

Council of Scientific and Industrial Research Fellowship 1974-1976. 

Indian National Merit Scholarship 1972-1974. 

Thatte memorial gold medal in Bachelor's degree 1972. 

PROFESSIONAL ACTIVITIES (National Committees and Forums): (Partial List, 1997-03 
only) 

1. Member, National Academies of Sciences Future of Supercomputing Committee, 2003- 
2004. 

2. Volume Editor, Wiley Encyclopedia of Genomics, Genetics, Proteomics and 
Bioinformatics. 2004. 

3. Member, Scientific Advisory Panel, NURSA project, NIDDK, NIH. 2004-present. 

4. Member, NIH Director's Roadmap Committee for Bioinformatics, 2003 

5. Chair, Software Maintenance and Development Study Section, NIH, 2002- present 

6. Chair, San Diego Supercomputer Center Executive Committee, 2002 - present 

7. Guest Professor, Centre for Molecular Biology and Neuroscience, University of Oslo, 
Norway 2002 - present 

8. Editorial Board Member, Bioinformatics. 2002-present. 

9. Editorial Board Member. Molecular and Cellular Proteomics 2002-present. 

10. Editorial Board Member. Journal of Molecular Recognition. 2000-present. 

11. Editor, Wiley Encyclopedia on Genomics, Proteomics and Bioinformatics 2002 - present 

12. Member, Search Committee for Director of Bioinformatics and Computational Biology 
Institute, NIGMS, NIH, 2001 

13. Chair, BISTI Study Section, NIH 2001-2002 

14. Chair, NIH/NSF Bioengineering and Bioinformatics Education and Training Workshop 
2001 

15. US Representative in the Global Science Forum Neuroinformatics Working Group 2000- 
present 

16. Member, Scientific Advisory Board to the Department of Bioengineering, Ul Chicago 
2000-present 

17. Member, Scientific Advisory Board, Mitokor Corporation, 2001-present 

18. Member, Scientific Advisory Board, Genomar, 1999-present 

19. Member, SUN Microsystems Informatics Advisory Council 2000-present 

20. Member, Review Panel, NSERC, Berkeley 1999-present 

21. Member, Advisory committee to the Working Group of the Director of NIH, on 
Bioinformatics (Chairs: Botstein/Smarr) 1998-99 

22. Member, NIH Site Visit Panel on Program Project Grant, Lawrence Berkeley Laboratory, 
1999. 

23. Member, Whitaker Foundation Workshop on Bioengineering Education. 1999 

24. Invited Participant, AlChe Planning Group in Bioinformatics, 1999. 

25. Adhoc Member, Genetics Study Section, NIH, 1998 

26. Panelist, Advisory Panel on Postdoctoral Fellowships in Bioinformatics of the National 
Science Foundation, 1999-2000 

27. Panelist, Committee on Knowledge and Distributed Intelligence Initiative, NSF 1998 
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28. Member, NIH Review Panel on Workshop and Training Grants in Computational 
Genomics and Bioinformatics, 1999. 

29. Chair, Session on Computational Biophysics Session, Annual Meeting of the Biophysical 
Society, 1999. 

30. Proposal Reviewer, W.M. Keck Foundation, 1998 

31. Member, Scientific Advisory Board, BioSoft Inc. Norway 1998-present 

32. Member, Advisory Committee, Fralin Biotechnology Center, Virginia Tech. 1998-99 

33. Chair & Organizer, Beckman Institute Symposia on Protein Structure, Function and 
Bioinformatics, 1991-97. - " ' ^ 

34. Member, Editorial Board for Journal of Molecular Recognition, 1999-present 

35. Ad-Hoc Member, SBIR/STTR Study Section, NIH, 1997-98. 

36. Member of the Working Group of the Advisory Committee to the Director of NIH on 
Biological Information Science and Technology Initiative. Participated in writing the 
BISTI report to NIH (1997-1999) 



INVITED LECTURES (Partial List 1991-1999): 

1. International Symposium on Polymer Modeling, NCSA. "Computer Simulations of 
Biomolecules". May 1991. 

2. Department of Chemistry, Northern Illinois University, DeKalb. "Molecular Recognition 
in Proteins". Oct. 1991. 

3. Eli Lilly Laboratories, Indianopolis. "Protein Structure Modeling". Sept. 1991. 

4. Iowa Academy of Sciences, Cedar Falls. "Supercomputing in Biology". Nov. 1991. 

5. Fourth International meeting on Software Engineering and Knowledge Engineering, 
Capri, Italy . "Knowledge-Based Approaches to Protein Structure and Motifs". June 
1992. 

6. Neural Networks - Biology to High Energy Physics, Elba, Italy. "Machine Learning 
Approaches to Protein Feature Prediction". June 1992. 

7. Department of Chemical Engineering, Illinois. "Computer Simulations of Biomolecular 
Recognition". Sept. 1992. 

8. Pfizer Central Research Laboratories. "Biomolecular Recognition", and " Protein 
Structure Prediction from Sequence". Sept. 1993. 

9. Center for Advanced Research in Biotechnology, University of Maryland. "Electrostatics 
and Molecular Recognition". Oct. 1993. 

10. International Symposium on Distance-based Approaches to Protein Structure, Lyngby, 
Denmark. Plenary Lecture on "Protein Structure Prediction - Past and Present". Nov. 
1993. 

11. North Dakota State University, Fargo. "Protein Structure and Modeling". April 1994. 
12. 1994 First World Congress on Computational Medicine, Health and Biotechnology. 

Session Chair and Speaker. "Knowledge-Based Approaches to Protein Structure". 
April 1994. 

13. Pfizer-Beckman Symposium on Protein Structure and Engineering, University of Illinois. 
"Knowledge-Based Tools for Protein Engineering and Design". June 1994. 

14. Antibody Workshop, Cambridge, U.K. "Kinetics and Energetics in Immune 
Recognition". Aug. 1994. 

15. Short Course on "Protein Modeling" in SISSA, Trieste, Italy. Sept. 1994. 

16. Texas A & M University, College Station. "Knowledge-Based Approaches to Protein 
Structure", Dec 1994. 
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17. University of California at Irvine. "Molecular Recognition in Proteins" - From Lock and 
Key to Float, Fix, Fit and Fasten., May 1995. 

18. ISIS Pharmaceuticals, San Diego, CA. "Protein Structure Prediction - Novel Methods"., 
May 1995. 

19. Intelligent Systems for Molecular Biology-95. Tutorial titled "Molecular Biology for the 
Computer Scientist". Cambridge, U.K. July 1995. 

20. Brandeis University, Dept. of Chemistry. "Knowledge-based Potentials for Proteins". 
Sept. 1995. 

21. Istituto Ricerca BioMoleculare P. Angeletti, Pomezia. "Protein Structure Prediction 
Methods". Nov. 1995. 

22. Centro Nazionale Ricerca (CNR) Genoa. "Simulations of Biological Membranes". Nov. 
1995. 

23. Rush Medical School, Chicago, IL. "Protein Structure and Folding". Jan. 1996. 

24. Univeristy of California at Irvine, "Protein Folding in machino". March. 1996. 

25. Purdue University, IN., "Knowledge-Based Methods for Protein Structure". March 1996. 

26. University of Wisconsin at Madison, Wl. "Molecular Recognition - Through a 
Computational Microscope". March 1996. 

27. International Center for Theoretical Physics, Trieste, Italy. "Methods for Protein 
Structure and Folding", March 1996. 

28. International Center for Theoretical Physics, Trieste, Italy. "Computer Simulation 
Methods for Molecular Recognition", March 1996. 

29. Eli Lilly and Company, Indianapolis. "The Biology Workbench ". April, 1996. 

30. ISMB96, Talk on "Knowledge-based Methods for Protein Structure Refinement and 
Prediction". June 1996. 

31. ISMB96, Presentation and National Release of the "Biology Workbench". June 1996. 

32. University of Minnesota, MN., "Protein Folding in machino" . Oct., 1996. 

33. University of Minnesota, MN., "Biology Workbench - A Seamless Integration of 
Dabatases, Analysis Algorithms and Interfaces on the World Wide Web". Oct. 1996. 

34. ASBMB Conference on Computational Biology - Methods in Biomolecular Imaging. "A 
Knowledge-Based Method for Protein Structure Refinement". Oct., 1996. 

35. USDA International Symposium on "Genetic Analysis of Economically Important Traits 
in Livestock" , Allerton, Illinois. "Biolnformatics on the World Wide Web". Nov., 1996. 

36. Chemical Engineering Symposium on Engineering Protein Recognition, Univ. of Illinois 
at Urbana-Champaign. "Role of Electrostatic Interactions in Antigen-Antibody 
Interactions". November, 1996. 

37. Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, 
NC. Workshop on Protein Folding. December, 1996. 

38. Midwest Regional ACS Meeting., Midland, ML, Bioinformatics on the Web., May 1997. 

39. NSF Workshop on "Knowledge Network Processing, U. Pennsylvania., June 1997. 

40. Second International Conference on Molecular Recognition on Surfaces, Basel, 
Switzerland, July 1997. 

41. Department of Chemical Engineering, University of Wisconsin at Madison, 
Bioinformatics and the Future., Sept. 1997. 

42. 55. Departments of Physics & Biochemistry, Oklahoma State University, Stillwater. 
Protein Folding in machino, Oct. 1997. 

43. Genome Workshop on Databases, "Bioinformatics and Genomics", Minnesuing Acres, 
Wl, Oct. 1997. 

44. Department of Biophysics, Johns Hopkins University, Protein Folding in machino, Dec. 
1997. 

45. National Institutes of Health BECON Symposium Panel on "Bioinformatics", Feb. 1997. 
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46. Institute for Biomedical Computing, Washington University, St. Louis, Protein Structure 
and Folding, March 1998. 

47. Departments of Computer Science and Biology, University of Kentucky, The Biology 
Workbench - A paradigm for virtual computing on the web, April 1998. 

48. Department of Biochemistry, Northern Illinois University, DeKalb. Molecular Recognition 
in the Immune System, April 1998. 

49. Department of Chemistry, Rice University. Protein Structure and Folding, April 1998. 

50. Department of Bioengineering, UC San Diego, La Jo||a. Bioinformatics - The Rosetta 
Stone of Biology, June 1998. 

51. Department of Bioengineering, UC San Diego, La Jolla. Molecular Recognition in 
Proteins, August 1998. 

52. European Molecular Biology Organization, "Talking Proteins Symposium", Heidelberg, 
Sept., 1998. 

53. Department of Computer Science, Purdue University, Bioinformatics., Oct., 1998. 

54. Department of Bioengineering, Rice University, Houston, Bioinformatics in the 
Postgenome Era, Nov., 1998. 

55. Chem. Eng. Dept. Massachusetts Institute of Technology, Bioinformatics: Deciphering 
the Rosetta Stone of Biology, April, 1999 

56. Chemistry Dept., Wabash College, Bioinformatics, April 1999 

57. Physiology & Biophysics, Cornell Medical School, Applications of Bioinformatics, NY, 
May 1999 

SINCE JOINING UCSD 

58. Symposium on "Bioinformatics, Patents and Intellectual Property Rights" organized by 
McDermott and Campbell & Flores Law Firms, La Jolla, San Diego, June 1999. 

59. Annual Symposium on Biomedical Engineering, "Bioinformatics in the Post-genome 
Era", Salt Lake City, Utah, July 1999. 

60. Chattaqua Conference on Computational Science, "Bioinformatics", Boston University, 
Sept. 1999. 

61. Department of Chemistry at UNC Chapel Hill, "Deconstructing Protein Structure", Oct., 
1999. 

62. Department of Chemistry and Biochemistry, Duke University, "Bioinformatics - 
Deciphering the Rosetta Stone of Biology", Oct., 1999. 

63. BMES Annual Symposium. Atlanta, GA, "Postgenomic Bioinformatics", Oct. 1999. 

64. Chair, Panel Discussion on Bioinformatics, BMES Annual Meeting, Atlanta, GA, Oct. 
1999. 

65. 3 rd TIGR Symposium on Computational Genomics, Baltimore, MD., "Sequence-Function 
Mapping - Lessons from Bioinformatics" , Nov., 1999. 

66. Department of Computer Science, San Diego State University, "Postgenomic 
Bioinformatics - Computer Science Challenges", Nov., 1999. 

67. CSU Biotechnology Symposium, Keynote Lecture, "Bioinformatics: Deciphering the 
Rosetta Stone of Biology", CalPolytechnic Pomona, Jan., 2000. 

68. BioStar Meeting on Bioinformatics, Bioinformatics: challenges and progress, UC San 
Diego., Jan 2000. 

69. Pacific Rim Society Internation Symposium on Bioinformatics. Bioinformatics at the 
Crossroads of Biology, Medicine and Information Science. Los Angeles, March, 2000. 

70. 19 th Southern Biomedical Conference, Keynote Lecture. The future of Bioinformatics., 

Blacksburg, Virginia, April 2000. 
71. Genetic Resources for the New Century, Beyond the Genome: Bioinformatics. San 

Diego Zoo., May 2000. 
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72. Global Forum on Neuroinformatics organized by National Institutes of Mental Health, 
"Design of a Neuroinformatics Workbench", Genoa, Italy, May 2000. 

73. Massachusetts Institute of Technology Workshop on Bioinformatics. "Protein Structure 
and Folding (4 Lectures)., MIT, Boston, June 2000. 

74. NIH Site Visit to evaluate a Program Project Grant, New Jersey Institute of Technology, 
June 2000. 

75. Iceland Genome Project. "Bioinformatics after the Genome" Reykjavik, July, 2000. 

76. International Society for Animal Genetics Conference, Keynote Lecture: "Bioinformatics: 
Mapping the New Central Dogma in Biology", Minneapolis, July 2000. 

77. Keynote Speaker, "The cross-education across biology and information science", 
Conference in Bioinformatics Education in Undergraduate Curricula. San Diego State 
University. July 2000. 

78. Mitokor Corporation, Designing Infrastructure for Bioinformatics. San Diego, August 
2000. 

79. Keynote Speaker, EUCHEM conference on BCDC- Bioinformatics, Cheminformatics, 
Datamining, MQSAR and Chemometrics, Stockholm, Sweden. Bioinformatics. Swedish 
Chemical Society, Stockholm. Septermber 2000. 

80. Mt. Sinai Medical School, Department of Pharmacology. Bioinformatics: Applications to 
sequence and function analysis. New York, October 2000. 

81. Conference, Metabolic Engineering III, Sequence-structure-function mapping for 
engineering microorganisms. Colorado Springs, October 2000. 

82. SUN Microsystems. Bioinformatics: Challenges for Software, Database and Computer 
Designers. Palo Alto, November 2000. 

83. PacificChem. Chair and Speaker, Bioinformatics Session. Beyond the Genome. 
December 2000. 

84. 5th Lake Tahoe Symposium on Molecular Diversity. Introduction to the Alliance for 
Cellular Signaling. Lake Tahoe, January 2001. 

85. Keynote and Award Lecture. LabAutomation 2001. Bioinformatics: At the Crossroads 
of Biology, Medicine and Engineering. Palm Springs. January 2001. 

86. UCSD Chemistry-Industry Presentation. Challenges in Post-Genome Bioinformatics. 
UCSD. February 2001. 

87. Health Care in the New Millenium: Merging Biology, Information Technology and 
Engineering, organized by AIMBE. The Information Revolution in Biology and its Impact 
on Medicine in the New Millenium. National Academy of Sciences. Washington D.C. 
March 2001. 

88. California State University Pomona, Departments of Chemistry and Biochemistry. 
Bioinformatics: Flow of Information in Biology. March 2001. 

89. BMES-FASEB Meeting. Chair and Speaker. Bioinformatics Session. Orlando April 
2001. 

90. Institute for Theoretical Physics, UC Santa Barbara. Modeling biological processes. 
April 2001. 

91. Department of Bioengineering. Ul Chicago. Advisory Committee Meeting. April 2001. 

92. Plenary Lecture. Symposium of European Society for Engineering and Medicine. 
Bioinformatics: The Interface between Biology, Engineering and Medicine. Belfast, 
Ireland. May 2001. 

93. Bioinformatics Course in MIT Workshop in Bioinformatics. MIT, Boston. June 2001. 

94. Databasing the Brain. University of Oslo and European Commission and Global 
Science Forum. Chair and Speaker. Bioinformatics and Neuroinformatics: Challenges 
and Opportunities. Oslo, Norway. July 2001. 

95. Salk Institute. Bioinformatics of Cellular Signaling. La Jolla. July 2001 . 
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96. Gordon Research Conference on Bioinformatics- From Inference to Predictive Models. 
Designing Knowledge Bases for Cellular Signaling. Tilton School, New Hampshire, 
August 2001. 

97. Aspen Workshop on Enzymatic Networks and Cellular Signaling. Data Modeling and 
Ontologies for AfCS. Aspen, August 2001 . 

98. Computational Challenges in the Post-Genomic Age II. Databases for Modern Biology. 
Keynote Lecture. Research Triangle Park, North Carolina. September, 2001. 

99. The O'Reilly Bioinformatics Technology Conference. State of the Art in Integrative 
Bioinformatics. Section Speaker. Tucson, Arizona. January 2002. 

100. Biophysical Society Annual Meeting. Semiotics of Cellular Signaling. San Francisco, 
February 2002 

101. University of Minnesota Bioinformatics Symposium. Challenges in Integrative 
Bioinformatics. Minneapolis, Minnesota, April 2002. 

102. In Silico Biology Conference. Molecular Ontology for Cellular Signaling. San Diego, 
June 2002. 

103. University of Utah Department of Medical Informatics. Postgenome Bioinformatics: At 
the Crossroads of Biology, Engineering and Medicine. September 2002. 

104. SIAM. Structuring Cellular Data for Computational Modeling. Washington DC, October 
2002. 

105. La Jolla Mesa Signaling Group. AFCS Molecule Pages. La Jolla, October 2002 

106. Oracle Life Sciences Conference. Structuring Biological Data. San Francisco, November 
2002 

107. International Business Forum. Investing and Partnering in Systems Biology. San 
Francisco, November 2002 

108. Harvard-MIT Seminar in Computational Biology. Challenges in Post Sequence 
Bioinformatics and Systems Biology. Boston, November 2002. 

109. Mount Sinai Scho.ol of Medicine. Modeling Cellular Complexity. New York City, 
December 2002. 

110. New York Academy of Sciences Conference. Understanding the Cell: Components, 
Networks and Modeling. New York City, December 2002. 

111. IBM Academy of Technology. Challenges in Post-Sequence Bioinformatics and 
Systems Biology. San Jose, California, March 2003 

112. Knowledge Millennium III: The Business of Biotechnology. New Delhi, India 

113. Purdue University. Systems Biology: Windows into Intracellular Networks. Lafayette, 
April 2003. 

114. Purdue University. Systems Biology: Challenges in the Post Genome Era: The Biology- 
Computational Science Interface. Lafayette, April 2003. 

1 15. * Virginia Bioinformatics Institute. Systems Biology: Challenges in the Post-Genomic Era. 

Bethesda, May 2003 

116. NIH Human Brain Project. Postgenomic Biology: Deciphering Intracellular Networks. 
Bethesda, May 2003 

117. Salk Institute. Deciphering Intracellular Networks in Mammalian Cells. La Jolla, August 
2003. 

118. Gordon Research Conference. Reconstruction and Modeling of Networks from Cellular 
Data. Queen's College, Oxford, August 2003. 

119. Johns Hopkins University, CCBM Inaugural Seminar, "Deciphering Intracellular 
Networks in Mammalian Cells". Baltimore, September 2003. 

120. CallT 2 . Systems Biology, Bioengineering and Bioinformatics: A new discipline at the 
Interface of Biology, Medicine and Information Technology. La Jolla, September 2003. 

121. Pathways, Networks and Systems Biology Conference. Deciphering Intracellular 
Networks in Mammalian Cells. Santorini, Greece, September 2003. 



Shankar Subramaniam 



Page 10 



122. Scientific Data Integration: Challenges and Some Solutions. Biomedical Information 
Science and Technology Initiative Consortium. Bethesda, MD, November 2003. 

123. 2003 PKR Protein Phosphorylation Workshop. .Alliance for Cellular Signaling: Design of 
the Informatics Infrastructure. Asilomar, CA December 2003. 

124. IPAM: Proteomics Workshop I: High Throughput Technologies and Methods of Analysis. 
"Reconstruction of Cellular Networks: the Biological Context". Los Angeles, March 2004. 

125. NHLBI: A Systems Biology Approach to Regulatory Networks in Heart, Lung, Blood and 
Sleep Research Working Group. " Biological .DataJntegration = and Analysis: Challenges 
for Mammalian Cell Systems. Bethesda, MD, May 2004. 

126. Massachusetts Institute of Technology: MIT Summer Course, "Bioinformatics: 
Principles, Methods and Applications" (five lectures). Boston MA, June 2004 

127. Stazione Zoologica: "Deciphering Mammalian Intracellular Networks: Integrative Data 
Analysis". Naples, Italy, July 2004. 

128. Iowa State University: "Cellular Systems Biology: Data Integration, Analysis and 
Modeling". Ames IA, July 2004 

129. University of Oslo: "Ontologies, Databases and Tools for Brain Imaging. Research 
Course. Oslo Norway. August 2004. 

130. IEEE-EMBC 04 Plenary Lecture: "Bioinformatics and Computational Systems Biology: 
At the Cross Roads of Biology, Engineering and Computation". San Francisco CA, 
September 2004. 

131. Science & Technology Expert Partnership (STEP) program, Conference on 
Computational Biology: "Data Integration and Analysis for the Reconstruction of 
Mammalian Cellular Networks". Washington DC, September 2004. 

132. Plenary Lecture, Frontiers in Oncology and Pathology Informatics, "Decoding networks 
in mammalian cells", Pittsburgh, October 2004. 

133. Special Lecture at the Nanotechnology Center, Purdue University, "Bioinformatics and 
Systems Biology: At the interface of biology, engineering and computer science". 
Purdue, October 2004. 

PUBLICATIONS: 

Refereed Publications 2001-2004 

1. Avidor-Reiss T, Maer AM, Koundakjian E, Polyanovsky A, Keil T, Subramaniam S, 
Zuker CS. Decoding cilia function; defining specialized genes required for 
compartmentalized cilia biogenesis. Cell. 117(4):527-39, 2004. 

2. Bornheimer, S.J., Maurya, M.R., Farquhar, M.G., and Subramaniam, S. Computational 
modeling reveals how interplay between components of a GTPase-cycle module 
regulates signal transduction. Proc. Natl. Acad. Sci. USA 101: 15899-15904, 2004. 

3. Hsaio A, Worrall DS, Olefsky JM and Subramaniam S. Variance-modeled posterior 
inference of microarray data: Detecting gene-expression changes in 3T3-L1 adipocytes. 
Bioinformatics. 2004 Nov 22:20(1 7):31 08-27. 

4. Papin, J. and Subramaniam, S. Bioinformatics and Cellular Signaling. Curr Opin 
Biotechnol. 1:78-81,2004. 

5. Chittababu Guda, Eoin Fahy and Shankar Subramaniam. MITOPRED: a genome-scale 
method for prediction of nucleus-encoded mitochondrial proteins. Bioinformatics 20:1-10 
2004. 

6. Giridhar Chukkapalli, Chittibabu Guda and Shankar Subramaniam. SledgeHMMER: a 
web server for batch searching the Pfam database. Nucleic Acids Research 32, Web 
Server issue 2004. 
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7. Chittibabu Guda, Purnima Guda, Eoin Fahy and Shankar Subramaniam. MITOPRED: a 
web server for the prediction of mitochondrial proteins. Nucleic Acids Research 32 2004. 

8. Jason Papin and Shankar Subramaniam. Bioinformatics and cellular signaling. Current 
Opinion in Biotechnology 15:78-81. 2004 

9. Dawn Cotter, Purnima Guda, Eoin Fahy and Shankar Subramaniam. MitoProteome: 
Mitochondrial Protein Sequence Database and Annotation System. Nucleic Acids 
Research 32: D463-D467. 2004. 

10. Robin T. Shealy, Anuradha D. Murphy, Rampriya Ramarathnam^Eric Jakobsson and 
Shankar Subramaniam. Sequence-Function Analysis of the K+-Selective Family of Ton 
Channels Using a Comprehensive Alignment and the KcsA Channel Structure. Biophys 
J 84:2929-2942, 2003. 

11. Dennis R. Livesay, Per Jambeck, Atipat Rojnuckarin and Shankar Subramaniam. 
Conservation of Electrostatic Properties within Enzyme Families and Superfamilies. 
Biochemistry 42:3464-3473 2003. 

12. Cornelius G. Hunter and Shankar Subramaniam. Protein Fragment Clustering and 
Canonical Local Shapes. Proteins: Structure, Function and Genetics. 50:580-588 2003. 

13. Cornelius G. Hunter and Shankar Subramaniam. Protein Local Structure Prediction from 
Sequence. Proteins: Structure, Function and Genetics. 50:572-579 2003. 

14. Subramaniam S, 2002, Bioinformatics of Cellular Signaling, Bock G and Goode J, eds, 
7n silico' simulation of Biological Processes, Wiley, Chichester (Novartis Foundation 
symposium 247), p 104-118. 

15. Joshua Li, Yuhong Ning, Warren Hedley, Brian Saunders, Yongsheng Chen, Nicole 
Tindill, Timo Hannay and Shankar Subramaniam. The Molecule Pages database. 
Nature 420:716-717. 2002. 

16. Gilman AG, Simon ML Bourne HR, Harris BA, Long R. Ross EM, Stull JT, Taussig R, 
Bourne HR, Arkin AP, Cobb MH, Cvster JG, Devreotes PN, Ferrell JE, Fruman D, Gold 
M, Weiss A, Stull JT, Berridge MJ, Cantlev LC, Catterall WA, Coughlin SR. Olson EN, 
Smith TF, Brugge JS, Botstein D, Dixon JE, Hunter T, Lefkowitz RJ, Pawson AJ, 
Sternberg PW, Varmus H, Subramaniam S, Sinkovits RS, Li J, Mock D, Ning Y, 
Saunders B, Sternweis PC, Hilgemann D, Scheuermann RH, DeCamp D, Hsueh R, Lin 
KM, Ni Y, Seaman WE, Simpson PC, O'Connell TP, Roach T, Simon Ml. Choi S, 
Eversole-Cire P, Fraser I. Mumbv MC, Zhao Y, Brekken D, Shu H, Meyer T, Chandv G, 
Heo WD, Liou J, O'Rourke N. Verghese M. Mumbv SM, Han H, Brown HA, Forrester JS, 
Ivanova P, Milne SB, Casey PJ, Harden TK, Arkin AP, Doyle J, Gray ML, Meyer T, 
Michnick S, Schmidt MA, Toner M, Tsien RY, Nataraian M, Ranganathan R, Sambrano 
GR; Participating investigators and scientists of the Alliance for Cellular Signaling. 
Overview of the Alliance for Cellular Signaling. Nature. 420(691 6):703-6, 2002. 

17. Eckerslev P, Egan GF, Amah S, Beltrame F. Bennett R, Biaalie JG, Dalkara T, De 
Schutter E, Gonzalez C, Grillner S, Herz A, Hoffmann KP, Jaaskelainen IP, Koslow SH, 
Lee SY, Matthiessen L, Miller PL, da Silva FM, Novak M, Ravindranath V, Ritz R. 
Ruotsalainen U, Subramaniam S, Toga AW, Usui S, van Pelt J, Verschure P, Willshaw 
D, Wrobel A, Tang Y; OECD Working Group on Neuroinformatics. Neuroscience data 
and tool sharing: a legal and policy framework for neuroinformatics. Neuroinformatics. 
1(2): 149-65,2003. 

18. Shun-lchi Amah, Francesco Beltrame, Jan G. Bjaalie, Turgay Dalkara, Erik De Schutter, 
Gary F. Egan, Nigel H. Goddard, Carmen Gonzalez Sten Grillner, Andreas Herz, K.- 
Peter Hoffmann, NroJaaskelainen, Stephen H. Koslow, Soo-Young Lee, Line 
Matthiessen, perry L. Miller, Fernanando Mira Da Silva, Mirko Novak, Viji Ravindranath, 
Raphael Ritz, Ulla Ruotsalainen, Vaclav Sebestra, Shankar Subramaniam, Yiyuan 
Tang, Arthur W. Toga, Shiro Usui, Jaap Van Pelt, Paul Verschure, David Willshaw and 
Andrzej Wrobel. Neuroinformatics: The Integration of Shared Databases and Tools 
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Towards Integrative Neuroscience. Journal of Integrative Neuroscience Vol. 1 No. 2 
(2002) 117-128. 

19. C. G. Hunter and S. Subramaniam. A Natural Coordinate Representation for the Protein 
Backbone Structure. Proteins: Structure, Function and Genetics. 49:206-215. 2003. 

20. R. J. Mashl, H. L. Scott, Subramaniam and E. Jakobsson. Molecular Simulation of 
Dioleoylphosphatidylcholine Lipid Bilayers at Differing Levels of Hydration. Biophys. J. 
81:3005-3015, 2001 

21. M. Farnum and S. Subramaniam. Assessment of CASP Models by a Knowledge-based 
Potential. Submitted to Proteins: Structure, Function and Genetics. 2001. 

22. Y. Fan and S. Subramaniam. Structure-specific Scoring Matrices. Submitted to J. 
Comput. Biology. 

23. D. R. Livesay and S. Subramaniam. Conserved sequence and structural association 
motifs in antibody-protein and antibody-hapten complexes. Submitted to Mol. 
Immunology. 2001. 

24. D. R. Livesay and S. Subramaniam. Electrostatics as the leit-motif of structural evolution 
in the CuZn superoxide dismutase protein family. Submitted to Biophys. J. 2001. 

25. M. Farnum, R. Leary and S. Subramaniam. Statistical Potentials and Folding Funnels 
for Protein Energy Landscapes. Manuscript in Preparation. 

26. J. Schnitzer and S. Subramaniam. Association of Cytochrome P450 and Ferrodoxins: 
Structure of Complexes from Brownian Simulations. Manuscript in Preparation. 

27. Y. Fan and S. Subramaniam. A Novel Strategy for Family-Specific Sequence Scoring 
Matrices. Manuscript in preparation. 

28. Bornheimer, S.J., Maurya, M., Farquhar, M. and Subramaniam, S. Computational 
Modeling of the GTPase Cycle and its Regulation. Submitted to Proceedings of National 
Academy of Sciences. 

29. Maurya, M., Bornheimer, S.J., Venkatasubramaniam, V., and Subramaniam, S. A 
multidimensional sensitivity analysis based framework for reduced order modeling: 
GTPase cycle as an exemplar. Submitted to Biophysical Journal. 

30. A. Hsiao, D.S. Worrall, J.M. Olefsky, and S. Subramaniam. Variance-modeled posterior 
inference of microarray data: Gene-expression changes in 3T3-L1 adipocytes Submitted 
to Bioinformatics. 

31. Ramarathnam, R. and Subramaniam, S. Phylogenomics of orthologous protein families: 
comparison of evolutionary profiles. Submitted for publication. 

32. Ramarathnam, R. and Subramaniam, S. Are Evolutionary Relationships in Proteins 
Driven by Complete Sequences or by Motifs? Submitted for Publication. 



Refereed Publications 1999-2000 

33. R. Ramarathnam and Subramaniam, S. A novel microarray strategy for detecting genes 
and pathways in microbes with unsequenced genomes. Microbial and Comparative 
Genomics. 5: 1-9, 2000. 

34. Herrgard S, Jambeck P, Taylor SS, Subramaniam S. Domain architecture of a 
caenorhabditis elegans AKAP suggests a novel AKAP function. FEBS Lett. 486:107-11, 
2000. 

35. G.Zou, R.D. Skeel, and S. Subramaniam. Biased Brownian Dynamics for Rate Constant 
Calculation. Biophys. J. 79: 638-645, 2000. 

36. Bouzat JL, McNeil LK, Robertson HM, Solter LF, Nixon JE, Beever JE, Gaskins HR, 
Olsen G, Subramaniam S, Sogin ML, Lewin HA. Phylogenomic Analysis of the Alpha 
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Proteasome Gene Family from Early Diverging Eukaryotes, J. Mol. Evol. 51:532-43, 
2000. 

37. A. Rojnuckarin, D. Livesay and S. Subramaniam. Bimolecular reaction simulation using 
weighted ensemble brownian dynamics and the university of houston brownian 
dynamics program. Biophys. J. 79: 686-693, 2000. 

38. S. Herrgard, C. Gibas and S. Subramaniam, Role of an electrostatic network of residues 
in the enzymatic action of Rhizomucor mihei lipase family. Biochemistry 39: 2921-2930 
2000. - - - - ^ - 

39. C. Gibas, P. Jambeck and S. Subramaniam. Role of pH in antibody-antigen 
complexation, Methods, Academic Press, Ed. (Linthicum & Subramaniam) 20: 292-309, 
2000. 

40. M. Viswanathan and S. Subramaniam. Molecular Dynamics Simulations of Antibody- 
Hapten Complexes, Methods, Academic Press, Ed. (Linthicum & Subramaniam) 20: 
362-371,2000. 

41. G. Altobelli & S. Subramaniam. Kinetics of Association of Anti-lysozyme Monoclonal 
Antibody D44.1 and Hen-egg Lysozyme. Biophys. J. .79:2954-65, 2000. 

42. S.W. Chiu, M. Clark, S. Subramaniam and E. Jakobsson. Collective Motion Artefacts 
Arising in Long-Duration Molecular Dynamics. J. Comp. Chem. 21:121-131, 2000. 

43. M. Viswanathan, D.S. Linthicum and S. Subramaniam. Analysis of Correlated Motion in 
Antibody Combining Sites from Molecular Dynamics Simulations, Methods 20: 362-371, 
2000. 

44. D. Livesay and S. Subramaniam. pH dependence of antibody-hapten association. 
Molecular Immunology 36: 397-410, 1999. 

45. Chiu SW, Jakobsson E, Subramaniam S, Scott HL. Combined monte carlo and 
molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl 
phosphatidylcholine lipid bilayers. Biophys. J. 77:2462-9, 1999. 

46. M.E. Wall, S. Subramaniam and G.N. Phillips. Protein structure determination using a 
database of interatomic distance probabilities. Protein Science 8:2720-2727, 1999. 

47. Rojnuckarin and Subramaniam, Knowledge-based potentials for protein structure, 
Proteins - Structure, Function and Genetics, 36: 54-67, 1999. 

48. S. W. Chiu, S. Subramaniam and E. Jakobsson, Simulation Study of Gramicidin/Lipid 
Bilayer System "in excess water and lipid I. Structure of the Molecular Complex, Biophys. 
Journal, 76:1929-38, 1999. 

49. S. W. Chiu, S. Subramaniam and E. Jakobsson, Simulation Study of Gramicidin/Lipid 
Bilayer System in excess water and lipid II. Rates and Mechanisms of Water Transport, 
Biophys. Journal, 76:1939-50, 1999. 

50. S.W. Chiu, M. Clark, E. Jakobsson, S. Subramaniam and H.L. Scott. Optimizations of 
Hydrocarbon Chain Interaction Parameters" Applications to the Simulation of Fluid 
Phase Lipid Bilayers. J. Phys. Chem. 103: 6323-6327, 1999. 

51.S.-W. Chiu, M. Clark, E. Jakobsson, S. Subramaniam, and H.L. Scott. Application of a 
Combined Monte Carlo and Molecular Dynamics Method to the Simulation of a 
Dipalmitoyl Phosphatidylcholine Lipid Bilayer. J. Comp. Chem. 20: 1153-1164, 1999. 



Refereed Publications Pre-1999 (Partial List) 

52. S. Subramaniam. The Biology Workbench - A Seamless Database and Analysis 
Environment for the Biologist, in "Bioinformatics", Proteins 32: 1-2, 1998. 

53. A. Rojnuckarin, S. Kim, S. Subrmaniam. Brownian Dynamics Simulations of Protein 
Folding: Access to Milliseconds Time Scale and Beyond. Proc. Natl. Acad. Sci. USA. 
95:4288-4292, 1998. 
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54. S. Sivasankar, S. Subramaniam and D. Leckband. Direct Measurement of the pH- 
Dependent Electrostatic Properties of a Protein Surface. Proc. Natl. Acad. Sci. USA 
95: 12961-12966.1998. 

55. J. Xiong, S. Subramaniam and Govindjee. A knowledge-based three dimensional model 
of the Photosystem II reaction center of Chlamydomonas reinhardtii. Photosynthesis 
Research 56: 229-254, 1998. 

56. Jie Liang, Pamidinghantam Sudhakar, Herbert Edelsbrunner, Ping Fu and Shankar 
Subramaniam. Analytical Shape Computing in Macromolecules. I. Calculation of 
Surface Areas and Volumes in Proteins. Proteins 33:1-17 1998. 

57. Jie Liang, Pamidinghantam Sudhakar, Herbert Edelsbrunner, Ping Fu and Shankar 
Subramaniam. Analytical Shape Computing in Macromolecules. II. Inaccessible Voids 
and Interfacial Regions in Proteins. Proteins 33:18-29 1998. 

58. H.L. Scott, E. Jakobsson, S. Subramaniam. Simulations of Lipid Membranes with 
Atomic Resolution. Comput. Phys. 12:328-333,1998. 

59. C. Gibas, S. Subramaniam, J.A. McCammon, B. Braden, R.J. Poljak. Electrostatic and 
Structural Factors in Antibody-Antigen Interactions. Biochemistry 36(50): 15599-15614, 
1997. 

60. J. Liang and S. Subramaniam. Computing Molecular Electrostatics A Boundary- 
Element Method for Protein Electrostatics. Biophys. J. 73: 1830-1841, 1997. 

61. C. Gibas and S. Subramaniam. Rules for Design of a Soluble Bacteriorhodopsin. Prot. 
Eng. 10: 1175-1190, 1997. 

62. G. Pappas and S. Subramaniam. An Evaluation of Protein Secondary Structure 
Prediction Algorithms. Tech. In Prot. Chem. VIII: 783-794, 1997. 

63. J. Xiong, S. Subramaniam and Govindjee. Modeling of D1/D2 proteins of the 
photosystem II reaction center of the cyanobacterium synechocystis sp. PCC 6803 and 
its implications to herbicide and bicarbonate binding. Prot. Sci. 5: 2054-2073, 1996. 

64. C. Singh, R. Sankararamakrishnan, S. Subramaniam, and E. Jakobsson. Solvation 
and Water Permeation of a Putative Model for the Pore Region of the Voltage-Gated 
Sodium Channel Biophys. J. 71:2276-2288,1996. 

65. M. Viswanathan, D.W. Pledger, S.Y. Tetin, D.S. Linthicum and S. Subramaniam. 
Modeling the Structure of the Combining Site of an Antisweet Taste Ligand Monoclonal 
Antibody NC 10.14. Biopolvmers 39: 395-406, 1996. 

66. S. Subramaniam, D.K. Tcheng and J.M. Fenton. Knowledge-Based Methods for Protein 
Structure Refinement and Prediction. In Proceedings of the Fourth International 
Conference on Intelligent Systems in Molecular Biology . St. Louis, 1996, Ed. David 
States et al., AAAI Press, California, pp. 218-229, 1996. 

67. C. Gibas and S. Subramaniam. Explicit Solvent Models in Protein pKa Calculations. 
Bioph. J 71: 138-147,1996. 

68. C. Mandal, M. M. Anchin, S. Subramaniam and D. S. Linthicum. ABGEN: A 
knowledge-Based Authomated Approach for Antibody Structure Modeling Nature - 
BioTechnology 14(3): 323-328,1996. 

69. L.Tang, T. G. Ebrey, and S.Subramaniam. Modeling the Structure of Visual Pigments. 
Isr. J. Chem. 35, 198-209, 1995. 

70. T. loerger, L. Rendell and S. Subramaniam. Change of representation to improve 
protein fold-class prediction. Machine Learning Journal 21:151-176, 1995. 

71. R. E. Kozack, M.J. D'Mello and S. Subramaniam. Computer Modeling of Electrostatic 
Steering and Orientational Effects in Antibody-Antigen Association. Biophys. J. 68(3): 
807-814, 1995. 

72. S. W. Chiu, M. Clark, V. Balaji, S. Subramaniam, H.L. Scott and E. Jakobsson. 
Incorporation of Surface Tension into Molecular Dynamics Simulation of an Interface: A 
Fluid Phase Lipid Bilayer Membrane. Biophys. J. 69(4): 1230-1 245, 1995. 
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73. M. Viswanathan, J. M. Anchin, P. R. Droupadi, C. Mandal, D.S. Linthicum, and S. 
Subramaniam. Structural Predictions of the Binding Site Architecture for Monoclonial 
Antibody NC6.8 using Ligand Binding, Spectroscopy and Computer-Aided Molecular 
Modeling Bioph. J. 69(3): 741-753, 1995. 

74. Jakobsson, E., Subramaniam, S., and H. L. Scott. In: Strategic issues in molecular 
dynamics simulations of membranes in: Membrane Structure and Dynamics , K. Merz 
and B. Roux, eds. Birkhauser, Boston, 1995. 

75. M. Hoist, R.E. Kozack, F.Saied, and S. Subramaniam.. Treatment of Electrostatic 
Effects in Proteins: Multigrid-Based- Newton Iterative Method for Solution of the Full 
Nonlinear Poisson-Boltzmann Equation. Proteins: Structure, Function and Genetics . 
18(3): 231-245, 1994. 

76. S. Subramaniam. Protein Structure Prediction-Past and Present. In: Protein Structure 
by Distance Analysis, Ed. H. Bohr & S. Brunak, IOS Press, 1994, pp 3-14. 

77. M. Hoist, R. E. Kozack, F. Saied, S. Subramaniam. Protein Electrostatics: Rapid 
Multigrid-Based Newton Algorithm for Solution of the Full Nonlinear Poisson-Boltzmann 
Equation. J. BiomoL Struct. Dvn. H(6): 1437-1445, 1994. 

78. K. E. Forsten, R. E. Kozack, D. A. Lauffenberger, and S. Subramaniam. Numerical 
Solution of the Nonlinear Poisson-Boltzmann Equation for a Membrane-Electrolyte 
System. J. Phvs. Chem. 98: 5580-5586, 1994. 

79. S. P. Slagle, R. E. Kozack and S. Subramaniam. Role of Electrostatics in Antibody- 
antigen Association: Anti-hen Egg Lysozyme/Lysozyme Complex (HyHEL-5/HEL) J. 
Biomol. Struct. Dvn . 12: 439-456, 1 994. 

80. R. Susnow, C. Schutt ,H. Rabitz and S. Subramaniam. Conformational Analysis of 
Dipeptides - A Sensitivity Analysis Approach. J. Comp. Chem. 15(9): 947-962, 1994. 

81. J. M. Anchin, C. Mandal, C. Culberson, S. Subramaniam and D. S. Linthicum. 
Computer-aided Modeling of the Binding Site Architecture for Eight Monclonal 
Antibodies that Bind a High Potency Guanidinium Sweetener. J. Mol. Graphics 12: 257- 
266, 1994. 

82. M. Reczko, H. Bohr, S.Subramaniam, S. Pamidinghantam and A. Hatzigeorgiou. Fold- 
Class Prediction by Neural Network. In: Protein Structure Prediction-Past and Present. 
In: Protein Structure by Distance Analysis, Ed. H. Bohr & S. Brunak, IOS Press, 1994, 
pp 277-286. 

83. R.E. Kozack and S. Subramaniam. Brownian Dynamics Simulations of Molecular 
Recognition in an Antibody-antigen System. Protein Science 2: 915-926, 1993. 

84. K. R. Rodgers, C. Su, S. Subramaniam and T. G. Spiro. Hemoglobin T-R 
Quaternary Structure and Dynamics from Simultaneous Monitoring of Tyrosine and 
Tryptophan UV Resonance Raman Signals. J. Am. Chem. Soc . 114 : 3697-3709, 1992. 

85. D. Bassblino, R. E. Bruccoleri & S. Subramaniam. Modeling the Antigen Combining 
Site of an Anti-dinitrophenyl Antibody. Protein Science 1 : 1465-1476, 1992. 

86. L. J. Nell, J. A. McCammon and S. Subramaniam. Anti-insulin Antibody. Structure and 
Conformation I. Molecular Modeling and Mechanics. Biopolvmers 32: 11-21 1992. 

87. V. Lounnas, B.M. Pettitt, L. Findsen & S. Subramaniam. A microscopic View of Protein 
Solvation. J. Phvs. Chem . 96: 7157-7159, 1992. 

88. S. Subramaniam, D. Tcheng, K. Hu, H. Raghavan & L. Rendell. Knowledge 
Engineering for Protein Structure and Motifs - Design of a Prototype System. IEEE 
Computer Society Journal (Proceedings of the Fourth International Conference on 
Software Engineering and Knowledge Engineering, Capri, Italy), 1992, pp. 420-434. 

89. D. K. Tcheng and S. Subramaniam. Machine Learning Approaches to Protein Feature 
Prediction. Int. J. Neural Systems. 3: 183-193, 1992. 
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90. S. W. Chiu, L. K. Nicholson, M. T. Brenneman, S. Subramaniam, Q. Teng, J. A. 
McCammon, T.A. Cross and E. Jakobsson. Molecular Dynamics Computations and 
Solid State NMR of the Gramicidin Cation Channel. Biophvs. J. 60: 974-978, 1991. 

91. S. W. Chiu, E. Jakobsson, S. Subramaniam and J. A. McCammon. Time Correlation 
Analysis of Water Transport through Gramicidin Ion Channel. Biophvs. J. 60: 273-285 
1991. 

92. P. H. Kussie, J. Anchin, S. Subramaniam, J. A. Glazel & D. S. Linthicum. Molecular 
modeling of a monoclonal antimorphine antibody. J. Immunolo gy 146 : 4248-4257 
1991. 

93. J. Anchin, S. Subramaniam and S. Linthicum. Structure and Characterization of an 
Anti-haloperidol Antibody. J. Mol. Recoq . 4: 7-15, 1991. 

94. J. Shen, C. F. Wong, S. Subramaniam and J. A. McCammon. Partial Electrostatic 
Charges for the Active Center of Cu, Zn Superoxide Dismutase. J. Comput. Chem . 3: 
346-350, 1990. 

95. S. W. Chiu, S. Subramaniam, E. Jakobsson and J. A. McCammon. Water and 
polypeptide conformations in gramicidin channel-A molecular dynamics study. Biophvs. 
J, 56: 253-261, 1989. 

96. S. Subramaniam, J.A. McCammon and R. J. Bacquet. Probing Molecular Recognition 
using Simulation Methods. In: The Immune Response to Structurally Defined Proteins, 
Ed. S.Smith-Gill and E.Sercarz, pp. 169-176, Adenine Press, 1989. 

97. C. F. Wong, J. Shen, S. Subramaniam, C. Zheng, and J. A. McCammon. Molecular 
Dynamics Simulation of Protein Hydration. J. Mol. Liquids . 41_: 193-206, 1989. 

98. J. Shen, S. Subramaniam, C. F. Wong, and J. A. McCammon. Superoxide dismutase- 
Fluctuations in the structure and solvationof the active site channel studied by molecular 
dynamics simulation. Biopolymers 28: 2085-2096,1989. 

99. W. J. Mortier, S. K. Ghosh and S. Subramaniam. Electronegativity equalization 
method and calculation of atomic charges in molecules. J. Amer. Chem. Soc. 108 : 
4315, 1986. 

100. S. Subramaniam, W. J. Mortier and S. K. Ghosh. Calculation of atomic 
charges in large molecules. Ann. N.Y. Acad. Sci . 482 : 82- 85, 1986. 

101. J. Hermans and S. Subramaniam. The free energy of xenon binding to 
myoglobin from molecular dynamics simulation. Isr. J. Chem. 27: 225-227,1986. 
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We begin our consideration of the generation of metabolic energy 
with glycolysis, a nearly universal pathway in biological systems. Gly- 
colysis is the sequence of reactions that converts glucose into pyruvate with the 
concomitant production of a relatively small amount of ATP. In aerobic organ- 
isms, glycolysis is the prelude to the citric acid cycle and the electron 
transport chain, which together harvest most of the energy contained in 
glucose. Under aerobic conditions, pyruvate enters mitochondria, where 
it is completely oxidized to C0 2 and H 2 0. If the supply of oxygen is 
insufficient, as in actively contracting muscle, pyruvate is converted into 
lactate. Under anaerobic conditions, yeast transforms pyruvate into etha- 
noL The formation of ethanol and lactate from glucose are examples of 
fermentations. 

The elucidation of glycolysis has a rich history. Indeed, the develop- 
ment of biochemistry and the delineation of this central pathway went 
hand in hand. A key discovery was made by Hans Buchner and Eduard 
Buchner in 1897, quite by accident. They were interested in manufactur- 
ing cell-free extracts of yeast for possible therapeutic use. These extracts 
had to be preserved without using antiseptics such as phenol, and so they 
decided to try sucrose, a commonly used preservative in kitchen chemis- 
try. They obtained a startling result: sucrose was rapidly fermented into 
alcohol by the yeast juice. The significance of this finding was immense. 
The Buchners demonstrated for the first time that fermentation could occur outside 
living cells. The accepted view of their day, asserted by Louis Pasteur in 
1860, was that fermentation is inextricably tied to living cells. The chance 



Glycolysis — 

Derived from the Greek 
stem glyk-, "sweet," and the 
word lysis, "dissolution." 



Fermentation — 

An ATP-generating process in 
which organic compounds act 
as both donors and acceptors of 
electrons. Fermentation can 
occur in the absence of 0 2 . 
Discovered by Pasteur, who 
described fermentation as "la 
vie sans lair" ("life without 
air"). 



Opening Image: Three-dimensional structure of phosphofructokinase, a key enzyme in 
glycolysis. Fructose 6-phosphate (red), a substrate, is next to A DP (blue), which occupies 
the ATP-binding site. The bacterial enzyme is allosterically activated by A DP (green) at 
o> distant site. Only one of the four subunits of the tetrameric enzy me is shown here. 
[Drawn from 4pjk.pdb. T. Schirmer and PR. Evans. Nature 343(1990): 140. The 
image was generated by Dr. Anthony Nicholls using GRASP, described in A. Nicholls, 
£ Sharp, and B. Honig. Proteins 1 1(l991):28l.] 




Light micrograph of yeast cells. 
[Courtesy of Dr. Randy Schekman.] 



Enzyme — 

A term coined by Friedrich 
Wilhelm Kuhne in 1878 to 
designate catalytically active 
substances that had previously 
been called ferments. Derived 
from the Greek words en, "in," 
and zyrne, "leaven." 



discovery of the Buchners refuted this vitalistic dogma and opened the 
door to modern biochemistry. Metabolism became chemistry. 

The next critical contribution was made by Arthur Harden and William 
Young in 1905. They added yeast juice to a solution of glucose and found 
that fermentation started almost immediately but soon ceased unless in- 
organic phosphate was added. They deduced that inorganic phosphate be- 
came linked to a sugar, and they proceeded to isolate a hexose diphosphate 
which was later shown to be fructose 1,6-bisphosphate, a key intermediate 
in glycolysis. Furthermore, Harden and Young discovered that yeast juice 
lost its activity if it was dialyzed or heated to 50°C. However, activity was 
restored when inactive dialyzed juice was mixed with inactive heated juice 
Thus, activity depended on the presence of two kinds of substances: a 
heat-labile, nondialyzable component (called zymase) and a heat-stable 
dialyzable fraction (called cozymase). We now know that "zymase" consiste 
of a number of enzymes, whereas "cozymase" consists of metal ions 
adenosine triphosphate (ATP), adenosine diphosphate (ADP), and coen- 
zymes such as nicotinamide adenine dinucleotide (NAD + ). 

Studies of muscle extracts then showed that many of the reactions of 
lactic fermentation were the same as those of alcoholic fermentation. This 
was an exciting discovery because it revealed an underlying unity in biochemistry 
The complete glycolytic pathway was elucidated by 1940, largely because 
of the pioneering contributions of Gustav Embden, Otto Meyerhof, Carl 
Neuberg, Jacob Parnas, Otto Warburg, Gerty Cori, and Carl Cori. Glycoly- 
sis is also known as the Embden-Meyerhqf pathway. 
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CH 3 — C— COO" Lactate 
OH 



Glucose 



Pyruvate ' o 
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Figure 19-1 

Some fates of glucose. 



CH 3 — CH 2 OH Ethanol 



AN OVERVIEW OF KEY STRUCTURES AND REACTIONS 

Learning the sequence of events in a metabolic pathway is easier with a 
firm grasp of the structures of the reactants and an understanding of the 
types of reactions taking place. Glycolytic intermediates have either six or 
three carbons. The six-carbon units are derivatives of glucose and fructose. 
The three-carbon units are derivatives of dihydroxyacetone, glyceraldehyde, glyc- 
erate, and pyruvate. 



q~ o o- 

-c-o-p-o- -c-o-p-o- 

I " II 

^ o o 



Es *er Anhydride drides. 



CH 2 OH V V V 

I III" 

H-C-OH H-C-OH C= 0 

CH 2° H CH 2 OH CH 2 OH CH 3 

Dihydroxyacetone Glyceraldehyde Glycerate Pyruvate 

All intermediates between glucose and pyruvate are phosphorylated. The 
phosphoryl groups in these compounds are linked as either esters or anhy- 
drides. J 



Now let us look at some of the kinds of reactions that occur in glycolysis 
and in many other metabolic pathways: 

1. Phosphoryl transfer. A phosphoryl group is transferred from ATP to a 
glycolytic intermediate, or from the intermediate to ADP, by a kinase. 
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R— OH + ATP 



Kinase 



R— 0— P— O" + ADP + 

I 

0~ 



2. Phosphoryl shift. A phosphoryl group is shifted from one oxygen 
atom to another within a molecule by a mutase. 



OH 



-O— P=0 
I 

0 



« JL . ^ " Mutase I 

R-C-CH 2 0— P— O" . R— C-CH 2 OH 



H 



H 



3. Isomerization. A ketose is converted into an aldose, 
an isomerase. 



or vice versa, by 



CH 2 OH 

I 

' _ Isomerase 
C=0 , 

I 

R 



Ketose 



C— H 
I 

H— C-OH 
I 

R 

Aldose 



4. Dehydration. A molecule of water is eliminated by a dehydratase. 



H — C — OH , DehYdrataSe 
I 

H 
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C— 
II 
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5. Aldol cleavage. A carbon-carbon bond is split in a reversal of an aldol 
condensation by an aldolase. 



R 
I 

C=0 



HO 



C-H 
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H— C-OH 
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R' 
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FORMATION OF FRUCTOSE 1,6-BISPHOSPHATE FROM 
GLUCOSE 

We now start our journey down the glycolytic pathway. The reactions in 
Ais pathway take place in the cell cytosol. The first stage, which is the 
conversion of glucose into fructose 1,6-bisphosphate, consists of three 
^eps: a phosphorylation, an isomerization, and a second phosphoryla- 
tion reaction. The strategy of these initial steps in glycolysis is to trap the substrate 




2 fxvr\ 



Electron micrograph of a liver cell. 
Glycolysis takes place in the cytosol. 
[Courtesy of Dr. Anne Hubbard.] 
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in the cell and form a compound that can be readily cleaved into phosphorylated 
three-carbon units. ATP is subsequently extracted from the three-carbon 
units. 



CH 2 OH 




2 -0 3 POH 2 C 



H OH 
Glucose 




CH 2 OP0 3 2 - 



H HO 
H\ l/OH 

OH H 

Fructose 1,6-bisphosphate 



Glucose enters most cells through specific transport proteins and has 
one principal fate: it is phosphorylated by ATP to form glucose 6-phosphate. The 
transfer of the phosphoryl group from ATP to the hydroxyl group on C-6 
of glucose is catalyzed by hexokinase. 




CH 2 OP0 3 2 - 



'H 

OH H/\ + ATP 
HO\L_l/OH 

H OH 
Glucose 



Hexokinase 




+ ADP + H + 



OH 

OH 

Glucose 6-phosphate 



Phosphoryl transfer is a basic reaction in biochemistry. As was dis- 
cussed earlier (p. 485), kinases are enzymes that catalyze the transfer of a 
phosphoryl group from ATP to an acceptor. Hexokinase, then, catalyzes 
the transfer of a phosphoryl group from ATP to a variety of six-carbon 
sugars (hexoses), such as glucose and mannose. Hexokinase, like all other 
kinases, requires Mg 2 + (or another divalent metal ion such as Mn 2+ ) for activity. 
The divalent metal ion forms a complex with ATP. The structures of two 
possible Mg 2+ -ATP complexes are shown in Figure 19-2. 



Figure 19-2 

Modes of binding Mg 2_l 



to ATP. 



Mg 2 ' 



-o 



. "? ?' 

Adenosine — 0 — P— 0 — P— 0— P — 0~ 



p 



?> 



Adenosine — 0 — P — Of \i 



The next step in glycolysis is the isomerization of glucose 6-phosphate to 
fructose 6-phosphate. The six-membered pyranose ring of glucose 6-phosphate 
is converted into the five-membered furanose ring of fructose 6-phosphate. 
Recall that the open-chain form of glucose has an aldehyde group on C-l, 
whereas the open-chain form of fructose has a keto group on C-2. The 
aldehyde on C-l reacts with the hydroxyl group on C-5 to form the pyra- 
nose ring, whereas the keto group on C-2 reacts with the C-5 hydroxyl to 
form the furanose ring. Thus, the isomerization of glucose 6-phosphate 
to fructose 6-phosphate is a conversion of an aldose into a ketose. 




H OH 

Glucose 6-phosphate 



Phosphoglucose 
isomerase 



2 -0 3 POH 9 C 



-^^0> N CH 2 0H 

\H H0>1 
H \[ I/ OH 

HO H 
Fructose 6-phosphate 



The open-chain representations of these sugars show the essence of this 
reaction. 
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V H 

I 

H— C— OH 

I 

HO— C— H 
I 

H— C— OH 
I 

H— C— OH 

CH 2 OP0 3 2 " 

Glucose 6-phosphate 

(An aldose) 



Phosphoglucose 
isomerase 



CH,OH 

I 

HO— C— H 

I 

H— C— OH 
I 

H— C— OH 

CH 2 OP0 3 2 " 

Fructose 6-phosphate 

(A ketose) 



A second phosphorylation reaction follows the isomerization step. Fruc- 
tose 6-phosphate is phosphorylated by ATP to fructose l f 6-bisphosphate. This com- 
pound was formerly known as fructose 1,6-diphosphate. Zfophosphate means 
two separate phosphate groups, whereas diphosphate (as in adenosine 
diphosphate) means two phosphate groups joined by an anhydride bond. 
Hence, the name fructose 1,6-bisphosphate should be used. 



2 ~0,POH 2 C 



CH 2 OH 



\ H HO / 

OH H 
Fructose 6-phosphate 



Phosphofructokihase 
+ ATP > 



2 "0 3 POH 2 C 



CH 2 OP0 3 2 " 



\ H Hoy 
H XI [/ OH 

OH H 
Fructose 1,6-bisphosphate 



+ ADP + H + 



This reaction is catalyzed by phosphofructokinase, an allosteric enzyme. The 
pace of glycolysis is critically dependent on the level of activity of this 
enzyme, which is allosterically controlled by ATP and several other metab- 
olites (p. 493). 



FORMATION OF GLYCERALDEHYDE 3-PHOSPHATE 
BY CLEAVAGE AND ISOMERIZATION 

The second stage of glycolysis consists of four steps, starting with the 
splitting of fructose 1,6-bisphosphate into glyceraldehyde 3 -phosphate and 
dihydroxyacetone phosphate. The remaining steps in glycolysis involve three- 
carbon units rather than six-carbon units. 



CH 2 OP0 3 2 " 
I 

c=o 

I 

HO— C— H 

I 

H— C— OH 
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H— C— OH 

CH 2 OP0 3 2 " 

Fructose 
1,6-bisphosphate 



Aldolase 



CH 2 OP0 3 
I 
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I 

HO— C— H 
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H 



Dihydroxyacetone 
phosphate 



v 

I 

+ H— C— OH 

CH 2 OP0 3 2 - 



Glyceraldehyde 
3-phosphate 



^This reaction is catalyzed by aldolase. This enzyme derives its name from 

nature of the reverse reaction, an aldol condensation. 
^( Glyceraldehyde 3-phosphate is on the direct pathway of glycolysis. Dihy- 
droxyacetone phosphate is not, but it can be readily converted into glycer- 



Aldol condensation — 
The combination of two 
carbonyl compounds (e.g., an 
aldehyde and a ketone) to form 
an aldol (a /3-hydroxy-carbonyl 
compound). 



T 

R 2 — C— R 3 + C 

A R 
Ketone Aldehyde 



to 

R2 — A — R3 

H— A— OH 

4. 
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Figure 19-3 

Triose phosphate isomerase consists of 
a central core of eight parallel ft 
strands (red) surrounded by eight a 
helices (green). Connecting regions 
are shown in yellow. This structural 
motif, called an aft barrel, is also 
found in one of the domains of pyru- 
vate kinase. [After a drawing kindly 
provided by Dr. Jane Richardson.] 
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1 ,3-Bisphosphoglycerate 
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aldehyde 3-phosphate. These compounds are isomers: dihydroxyacetone 
phosphate is a ketose, whereas glyceraldehyde 3-phosphate is an aldose. 
The isomerization of these three-carbon phosphorylated sugars is cata- 
lyzed by triose phosphate isomerase (Figure 19-3). This reaction is rapid and 
reversible. At equilibrium, 96% of the triose phosphate is dihydroxyace- 
tone phosphate. However, the reaction proceeds readily from dihydroxy- 
acetone phosphate to glyceraldehyde 3-phosphate because of efficient 
removal of this product by subsequent reactions. 



CH 2 0H 

I 

?-° 

CH 2 OP0 3 2 - 

Dihydroxyacetone phosphate 

(A ketose) 



Triose phosphate 
isomerase 
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I 

H— C— OH 
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Thus, two molecules of glyceraldehyde 3-phosphate are formed from one 
molecule of fructose 1,6-bisphosphate by the sequential action of aldolase 
and triose phosphate isomerase. The economy of metabolism is evident 
in this reaction sequence. The isomerase funnels dihydroxyacetone into 
the main glycolytic pathway — a separate set of reactions are not needed. 



ENERGY CONSERVATION: PHOSPHORYLATION IS COUPLED 
TO THE OXIDATION OF GLYCERALDEHYDE 3-PHOSPHATE 

The preceding steps in glycolysis have transformed one molecule of glu- 
cose into two molecules of glyceraldehyde 3-phosphate. No energy has yet 
been extracted. On the contrary, two molecules of ATP have been in- 
vested thus far. We come now to a series of steps that harvest some of the 
energy contained in glyceraldehyde 3-phosphate. The initial reaction in 
this sequence is the conversion of glyceraldehyde 3-phosphate into 1,3- 
bisphosphoglycerate (1,3-BPG), a reaction catalyzed by glyceraldehyde 3-phos- 
phate dehydrogenase. In the earlier literature, 1,3-BPG was known as 1,3- 
diphosphoglycerate (1,3-DPG). 

A high-potential phosphorylated compound is generated in this oxidation- 
reduction reaction. The aldehyde group at C-l is converted into an acyl 
phosphate, which is a mixed anhydride of phosphoric acid and a carboxylic 



acid The energy for the formation of this anhydride, which has a high 
phtphor y r g roVtransfer potential, comes from the oxidation of the 
aSe group. Note that C I in 1,3-BPG is at the oxidation level of a 
caAoSic acid NAD- (p. 450) is the electron acceptor in this oxidation^ 
The mechanism of this complex reaction, which couples oxidat.on and 
phosphorylation, will be discussed later in this chapter (p. 501). 

FORMATION OF ATP FROM 1,3-BISPHOSPHOGLYCERATE 

In the next step, the high phosphoryl transfer potential of 1,3-BPG is 
used to generate ATP. Indeed, this is the first ATP-generating reaction in 
Zco^si Phosphoglycerate kinase catalyzes the transfer of the phosphoryl 
g rouXm the afyl phosphate of 1,3-BPG to ADP. ATP and 3-phospho- 

^Z l ~ C ot the reactions catalyzed by glyceraldehyde 
3-phosphate dehydrogenase and phosphoglycerate kinase are 

1. Glyceraldehyde 3-phosphate, an aldehyde, is oxidized to 3-phospho- 
glycerate, a carboxylic acid. 

2. NAD 4 " is concomitandy reduced to NADH. 

3. ATP is formed from Pi and ADP. 

FORMATION OF PYRUVATE AND THE GENERATION 
OF A SECOND ATP 

In the last stage of glycolysis, 3-phosphoglycerate is converted into pyru- 
vate, and a second molecule of ATP is formed. 
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• The first reaction is a rearrangement. The position of the phosphoryl 
•group shifts in the conversion of 3-phosphoglycerate into 2-phosphogycerate, a 
•reaction catalyzed by phosphoglycerate mutase. In general, a mutase is an 
.enzyme that catalyzes the intramolecular shift of a chemical group, such 

^^o»on, an enol is formed by the dehydration of 
S-phosphoglycerate. Enolase catalyzes the formation of phosphoenolpyru- 
vate. This dehydration reaction markedly elevates the group transfer po- 
tential of the phosphoryl group. An enol phosphate has a high phosphory 
transfer potential, whereas the phosphate ester of an ordinary alcohol has 
a low one. The reasons for this difference will be discussed later (p. 504) . 
V In the last reaction, pyruvate is formed, and ATP is generated concomi- 
tantly. The virtually irreversible transfer of a phosphoryl group from 
phosphoenolpyruvate to ADP is catalyzed by pyruvate kinase. 
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ENERGY YIELD IN THE CONVERSION OF GLUCOSE 
INTO PYRUVATE 



The net reaction in the transformation of glucose into pyruvate is 

Glucose + 2 Pj + 2 ADP + 2 NAD + > 

2 pyruvate + 2 ATP + 2 NADH + 2 H + + 2 H 2 0 
Thus, two molecules of ATP are generated in the conversion of glucose into two 
molecules of pyruvate. A summary of the steps in which ATP is consumed or 
formed is given in Table 19-1. Recall that a pair of three-carbon units are 
formed from fructose 1,6-bisphosphate. The reactions of glycolysis are 
summarized in Figure 19-4 and Table 19-2. 



Figure 19-4 

The glycolytic pathway. 
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Table 19-1 

Consumption and generation of ATP in glycolysis 



Reaction 



Glucose > glucose 6-phosphate 

Fructose 6-phosphate > fructose 1,6-bisphosphate 

2 1,3-Btsphosphoglycerate > 2 3-phosphoglycerate 

2 Phosphoenolpyruvate > 2 pyruvate 



ATP change per 
glucose 



-1 
-1 
+2 
+2 

Net +2 
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Table 19-2 

Reactions of glycolysis 



Step 


Reaction 


Enzyme 


Tvoe* 


AG°' 


AG 


1 


Glucose + ATP » glucose 6-phosphate + ADP + H + 


Hexokinase 


a 


-4.0 


-8.0 


2 


Glucose 6-phosphate fructose 6-phosphate 


Phosphoglucose isomerase 


c 


+0.4 


-0.6 


3 


Fructose 6-phosphate + ATP > 

fructose 1,6-bisphosphate + ADP + H + 


Phosphofructokinase 


a 


-3.4 


-5.3 


4 


Fructose 1,6-bisphosphate 7 — " 

dihydroxyacetone phosphate + glyceraldehyde 

3-phosphate 


Aldolase 


e 


+5.7 


-0.3 


, 5 


Dihydroxyacetone phosphate t — " glyceraldehyde 

3-phosphate 


Triose phosphate isomerase 


c 


+ 1.8 


+0.6 


6 


Glyceraldehyde 3-phosphate + P, + NAD + 7 — * 

1,3-bisphosphoglycerate + NADH +H + 


Glyceraldehyde 3-phosphate 
dehydrogenase 


f 


+ 1.5 


-0.4 


. 7 


1,3-Bisphosphoglycerate + ADP 7 — - 

3-phosphoglycerate + ATP 


Phosphoglycerate kinase 


a 


-4.5 


+0.3 


' 8 


3-Phosphoglycerate 2-phosphoglycerate 


Phosphoglycerate mutase 


b 


+ 1.1 


+0.2 


■ ' 9 ■- 


2-Phosphoglycerate phosphoenolpyruvate +H 2 0 


Enolase 


d 


+0.4 


-0.8 


10 


Phosphoenolpyruvate + ADP + H + > pyruvate +ATP 


Pyruvate kinase 


a 


-7.5 


-4.0 



tfReaction type: (a) phosphoryl transfer; (b) phosphoryl shift; (c) isomerization; 
' (d) dehydration; (e) aldol cleavage; (f) phosphorylation coupled to oxidation. 

Note: AG°' and AG are expressed in kcal/mol. AG, the actual free-energy change, has 
been calculated from AG 0 ' and known concentrations of reactants under typical physiologic 
.conditions. Glycolysis can proceed only if the AG values of all reactions are negative. The 
small positive AG values of three of the above reactions indicate that the concentrations of 
metabolites in vivo in cells undergoing glycolysis are not precisely known. 
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|NTRY OF FRUCTOSE AND GALACTOSE 
!!yto GLYCOLYSIS 



Fructose 



6§t us consider how two other abundant sugars — fructose and galactose — 
&|n be funneled into the glycolytic pathway. Recall that the hydrolysis of 
jjljfrose (common table sugar) yields fructose and glucose, and that hy- 
ysis of lactose (milk sugar) gives galactose and glucose. Fructose itself 
^present in many foods (e.g., honey). A typical daily intake of fructose is 
&9 grams. Much of it is metabolized by the liver, using the fructose 1-phos- 
P^te pathway (Figure 19-5). The first step is the phosphorylation of fruc- 
^se to fructose 1-phosphate by fructokinase. Fructose 1-phosphate is then 
lljt into glyceraldehyde and dihydroxyacetone phosphate. This aldol 
avage is catalyzed by a specific fructose 1-phosphate aldolase. Glyceral- 
ii^yde is then phosphorylated to glyceraldehyde 3-phosphate by triose 
^ e so that it too cau enter glycolysis. 
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Fructokinase 



Fructose 1-phosphate 



Fructose 1-phosphate 
afdofase 



, , , Dihydroxyacetone 
Glyceraldehyde + phosphate 



Triose 
kinase 



r ATP 



ADP 



Glyceraldehyde 3-phosphate 
Figure 19-5 

Fructose enters the glycolytic pathway 
via the fructose 1-phosphate pathway; 



Alternatively, ymctoe can be phosphorylated to fructose 6-phosphate by hexoki- 
nose. However, the affinity of hexokinase for glucose is 20 times as high as 
it is for fructose. Little fructose 6-phosphate is formed in the liver because 
of the abundance of glucose relative to fructose in this organ. In contrast, 
adipose tissue has much more fructose than glucose. Hence, the forma- 
tion of fructose 6-phosphate is not competitively inhibited to an appreci- 
able extent, and most of the fructose in adipose tissue is metabolized 
through fructose 6-phosphate. 

Galactose is converted into glucose 6-phosphate in four steps. The first 
reaction in the galactose -glucose interconversion pathway is the phosphoryla- 
tion of galactose to galactose 1-phosphate by galactokinase. 

Galactose + ATP > galactose 1-phosphate + ADP + H + 

Galactose 1-phosphate then acquires a uridyl group from uridine diphos- 
phate glucose (UDP-glucose), an intermediate in the synthesis of glyco- 
side linkages (p. 586) . The products of this reaction, which is catalyzed by 
galactose 1-phosphate uridyl transferase, are UDP-galactose and glucose 
1-phosphate. 
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The galactose moiety of UDP is then epimerized to glucose. The configu- 
ration of the hydroxyl group at C-4 is inverted by UDP-galactose 4-epimer- 
ase. The sum of the reactions catalyzed by galactokinase, the transferase, 
and the epimerase is 

Galactose + ATP > glucose 1-phosphate + ADP + H + 

Note that UDP-glucose is not consumed in the conversion of galactose 
to glucose because it is regenerated from UDP-galactose by the epimer- 
ase. This reaction is reversible, and the product of the reverse direction is 
also important. The conversion of UDP-glucose into UDP-galactose is essential 
for the synthesis of galactosyl residues in complex polysaccharides and glycoproteins 
if the amount of galactose in the diet is inadequate to meet these needs. 

Finally, glucose 1-phosphate, formed from galactose, is isomerized to 
glucose 6-phosphate by phosphoglucomutase. We shall return to this reac- 
tion when we consider the synthesis and degradation of glycogen, which 
proceeds through glucose 1-phosphate (p. 584). 



GALACTOSE IS HIGHLY TOXIC IF THE TRANSFERASE 
IS MISSING 
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The absence of galactose 1-phosphate uridyl transferase causes galactose- 
mia, a severe disease that is inherited as an autosomal recessive trait. The 
metabolism of galactose in people who have this disease is blocked at 
galactose 1-phosphate. Afflicted infants fail to thrive. Vomiting or diar- 
rhea occurs when milk is consumed, and enlargement of the liver and 
jaundice are common. Furthermore, many galactosemics become men- 
tally retarded. The blood galactose level is markedly elevated, and galac- 
tose is found in the urine. The absence of the transferase in red blood 
cells is a definitive diagnostic criterion. 

Galactosemia is treated by the exclusion of galactose from the diet. A galactose- 
free diet leads to a striking regression of virtually all the clinical symp- 
toms, except for mental retardation, which may not be reversible. Contin- 
ued galactose intake may lead to death in some patients. The damage in 
galactosemia is caused by an accumulation of toxic substances, rather than by the 
absence of an essential compound. Patients are able to synthesize UDP- 
galactose from UDP-glucose because their epimerase activity is normal. 
One of the toxic substances is galactitol, which is formed by reduction of 
galactose. The presence of aldose reductase in the lens of the eye causes 
. galactitol to accumulate there, which leads to the entry of water and the 
development of cataracts. 



PHOSPHOFRUCTOKINASE IS THE KEY ENZYME 
IN THE CONTROL OF GLYCOLYSIS 

The glycolytic pathway has a dual role: it degrades glucose to generate 
ATP, and it provides building blocks for synthetic reactions, such as the 
-formation of long-chain fatty acids. The rate of conversion of glucose into 
pyruvate is regulated to meet these two major cellular needs. In metabolic 
'pathways, enzymes catalyzing essentially irreversible reactions are potential sites of 
control. In glycolysis, the reactions catalyzed by hexokinase, phosphofruc- 
tokinase, and pyruvate kinase are virtually irreversible; hence, they would 
be expected to have regulatory as well as catalytic roles. In fact, each of 
them serves as a control site. Their activities are regulated by the revers- 
ible binding of allosteric effectors or by covalent modification. Also, the 
amounts of these key enzymes are varied by transcriptional control to 
meet changing metabolic needs. Reversible allosteric control, regulation 
*by phosphorylation, and transcriptional control typically occur in times 
of milliseconds, seconds, and hours, respectively. 

Phosphofructokinase is the most important control element in the glycolytic path- 
way of mammals. The enzyme from liver (a 340-kd tetramer) is inhibited by 
high levels of ATP, which lower its affinity for fructose 6-phosphate. A high" 
concentration of ATP converts the hyperbolic binding curve of fructose 
6-phosphate into a sigmoidal one (Figure 19-6). This allosteric effect is 
elicited by the binding of ATP to a specific regulatory site that is distinct 
from the catalytic site. The inhibitory action of ATP is reversed by AMP, 
and so the activity of the enzyme increases when the ATP/ AMP ratio is lowered. In 
other words, glycolysis is stimulated as the energy charge falls. A second control 
feature comes into play when the pH drops appreciably. The inhibition of 
phosphofructokinase by H + prevents excessive formation of lactate 
(p. 497) arid a precipitous drop in blood pH (acidosis). 

Glycolysis also furnishes carbon skeletons for biosyntheses, and so phos- 
phofructokinase should also be regulated by a signal indicating whether 
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Figure 19-6 

Allosteric regulation of phosphofruc- 
tokinase. A high level of ATP inhibits 
the enzyme by decreasing its affinity 
for fructose 6-phosphate. AMP dimin- 
ishes and citrate enhances the inhibi- 
tory effect of ATP. 
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Fructose 2,6-bisphosphate 
(F-2,6-BP) 



^ ab r dam ° r SCarCC - IndCed ' P^P^tokinase is in- 
hibited by citrate, an early intermediate in the citric acid cycle (p 510) A 

high level of citrate means that biosynthetic precursors are abundant and 
so additional glucose should not be degraded for this purpose. Citrate 
inhibits phosphofructokinase by enhancing the inhibitory effect of ATP 
A new regulator of glycolysis was discovered in 1980 by Henri-Gerv 
Hers and Emile Van Schaftingen. They found fructose 2,6-bisphosphal 
a previously unknown metabolite, is a potent activator of phospLfruct^ 
kinase. Fructose 2,6-bisphosphate (F-2,6-BP) activates phosphofructokt 
nase in liver by increasing its affinity for fructose 6-phosphate and dimin- 
ishing the inhibitory effect of ATP (Figure 19-7). In essence, F-2 6-BpI 
an alhsteric activator that shifts the conformational equilibrium of this tetrameric 
enzyme from the T state to the R state. 
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Figure 19-7 



A REGULATED BIFUNCTIONAL ENZYME SYNTHESIZES 
AND DEGRADES FRUCTOSE 2,6-BISPHOSPHATE 

fDhos S n e h a ; 6 " biSPh ° Sphate ^ f ° rmed by thC Phosphorylation of fructose 
b-phosphate, ,n a reaction catalyzed by phosphofructokinase 2 (PFK2) a 
different enzyme from phosphofructokinase (PFK) (Figure 19-8) F-2 6- 
BP is hydrolyzed to fructose 6- P hosphate by a specific pLsphataL fruc- 

Zf£7 T 2 (FBPaSC2) - Strik -/ fi ^ing P is tLt b ZhWK2 
ZiZ r ™.P rM ™ a ^gle 55-kd polypeptide chain. This bifunctional 
enzyme contains an N-terminal regulatory domain, followed by a kinase 
domain and a phosphatase domain. It is interesting to note that PF^ re em 
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bles phosphofructokinase, whereas FBPase2 resembles phosphoglycerate Chapter 19 495 

mutase. The bifunctional enzyme probably arose by the fusion of genes 
encoding the kinase and phosphatase domains. 

Fructose 6-phosphate accelerates the synthesis of F-2,6-BP and inhibits 
its hydrolysis. Hence, an abundance of fructose 6-phosphate leads to a higher 
concentration ofF-2,6-BP, which in turn stimulates phosphofructokinase. Such a 
process is called feedforward stimulation. Furthermore, the activities of 



Fructose 6-phosphate 



Stimulated by 
fructose 6-phosphate 




Inhibited by 
fructose 6-phosphate 



Fructose 2,6-bisphosphate 



Figure 19-8 

Control of the synthesis and degrada- 
tion of fructose 2,6-bisphosphate. 
Fructose 6-phosphate accelerates the 
formation of F-2,6-BP and inhibits its 
hydrolysis. A low blood glucose level 
leads to a higher level of the phos- 
phorylated bifunctional enzyme and 
hence to a lower level of F-2.6-BP. 



PFK2 and FBPase2 are reciprocally controlled by phosphorylation of a single 
serine residue. When glucose is scarce, a rise in the blood level of the hor- 
mone glucagon triggers a cyclic AMP cascade leading to the phosphoryla- 
tion of this bifunctional enzyme. This covalent modification activates 
FBPase2 and inhibits PFK2, lowering the level of F-2,6-BP. Conversely, 
when glucose is abundant, the enzyme loses its attached phosphate 
group, which leads to a rise in the level of F-2,6-BP and the consequent 
acceleration of glycolysis. This coordinated control is facilitated by having 
the kinase and phosphatase domains on the same polypeptide chain as 
the regulatory domain. We shall return to this elegant switch when we 
consider the integration of carbohydrate metabolism (p. 766). 

HEXOKINASE AND PYRUVATE KINASE ALSO SET THE PACE 
OF GLYCOLYSIS 



'Hexokinase, the enzyme catalyzing the first step of glycolysis, is inhibited 
by glucose 6-phosphate. When phosphofructokinase is inactive, the con- 
centration of fructose 6-phosphate rises. In turn, the level of glucose 
6-phosphate rises because it is in equilibrium with fructose 6-phosphate. 
Hence, the inhibition of phosphofructokinase leads to the inhibition of hexokinase. 

bwever, liver possesses glucokinase, a specialized isoform of hexokinase 
that is not inhibited by glucose 6-phosphate. Glucokinase phosphorylates 
|lucose only when it is abundant because it has a much higher K M for 
|lucose than does hexokinase (5 mM, compared with 0.1 itim). The role of 
"glucokinase is to provide glucose 6-phosphate for the synthesis of glyco- 
gen, a storage form of glucose (p. 586). The high K M of glucokinase in 
fhe liver gives brain and muscle first call on glucose when its supply is 
fenited. 

|: Why i s phosphofructokinase rather than hexokinase the pacemaker of 
|tycolysis? The reason becomes evident on noting that glucose 6-phos- 
Rhate is not solely a glycolytic intermediate. Glucose 6-phosphate can also 
g? converted into glycogen or it can be oxidized by the pentose phos- 
jpte pathway (p. 559) to form NADPH. The first irreversible reaction 




Structure of a bacterial phosphofruc- 
tokinase. The four identical subunits 
are shown in different colors in a 
model depicting the surface of the 
enzyme. [Courtesy of Dr. Anthony 
Nicholls. Drawn from 4pfk.pdb. 
T Sch inner and P.R. Evans. Nature 
343 (1990): 140.] 
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unique to the glycolytic pathway, called the committed step, is the phos- 
phorylation of fructose 6-phosphate to fructose 1,6-bisphosphate. Thus, it 
is highly appropriate for phosphofructokinase to be the primary control 
site in glycolysis. In general, the enzyme catalyzing the committed step in a meta- 
bolic sequence is the most important control element in the pathway. 

Pyruvate kinase, the enzyme catalyzing the third irreversible step in 
glycolysis, controls the outflow from this pathway. This final step yields 
ATP and pyruvate, a central metabolic intermediate that can be oxidized 
further or used as a building block. Several forms of pyruvate kinase (a 
tetramer of 57-kd subunits) encoded by different genes are present in 
mammals: the L type predominates in liver, and the M type in muscle and 
brain. These variations on a common theme, called isoenzymes or isozymes, 
have essentially the same architectural plan and catalytic mechanism but 
differ in how they are regulated. Both the L and M forms bind phospho- 
enolpyruvate cooperatively. Fructose 1,6-bisphosphate, the product of 
the preceding irreversible step in glycolysis, activates pyruvate kinase to 
enable it to keep pace with the oncoming high flux of intermediates. ATP 
allosterically inhibits pyruvate kinase to slow glycolysis when the energy 
charge is high (Figure 19-9). 



Figure 19-9 

Control of the catalytic activity of 
pyruvate kinase. 
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Alanine (synthesized in one step from pyruvate, p. 630) also allosteri- 
cally inhibits pyruvate kinase, in this case to signal that building blocks are 
abundant The catalytic properties of the L form— but not of the M form- 
are also controlled by reversible phosphorylation. When the blood glu- 
cose level is low, glucagon triggers a cyclic AMP cascade that leads to the 
phosphorylation of pyruvate kinase, which diminishes its activity. A rise 
in the cytosolic calcium level induced by hormones such as vasopressin 
also leads to phosphorylation and inhibition of pyruvate kinase. These 
hormone-triggered phosphorylations, like that of the bifunctional enzyme control- 
ling the levels of fructose 2,6-bisphosphate, prevent the liver from consuming glu- 
cose when it is more urgently needed by brain and muscle (p. 575) . We see here a 
clear-cut example of how isoenzymes contribute to the metabolic diver- 
sity of different organs. 



DIVERSE FATES OF PYRUVATE: 

ETHANOL, LACTATE, OR ACETYL COENZYME A 

The sequence of reactions from glucose to pyruvate is similar in all organ- 
isms and in all kinds of cells. In contrast, the fate of pyruvate is variable. 
Three reactions of pyruvate are of prime importance: 

1. Ethanol \s formed from pyruvate in yeast and several other microor- 
ganisms. The first step is the decarboxylation of pyruvate. 



Pyruvate + H + > acetaldehyde + C0 2 Chapter 19 497 

This reaction is catalyzed by pyruvate decarboxylase. Thiamine pyrophos- GLYCOLYSIS 
phate, the coenzyme here, also participates in reactions catalyzed by de- 
carboxylases (p. 516) and other enzymes (p. 566). 
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The second step is the reduction of acetaldehyde to ethanol by NADH, 
in a reaction catalyzed by alcohol dehydrogenase. 

Acetaldehyde + NADH + H + ethanol + NAD+ 

; The active site of alcohol dehydrogenase contains a zinc ion that is coordi- 
nated to the sulfur atoms of two cysteine residues and a histidine nitrogen 
I: atom (Figure 19-10). As in carboxypeptidase A (p. 221), Zn 2+ polarizes 
; the carbonyl group of the substrate to stabilize the transition state (Figure 
19-11). 




Figure 19-10 

Model of the active site of liver alco- 
hol dehydrogenase showing the coor- 
dination of the zinc ion. [Drawn from 
2ohx.pdb; coordinates deposited by 
S. Al-Karadaghi and E.S. Cedergren- 
Zeppezauer.] 
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Figure 19-11 

Catalytic mechanism of alcohol dehy- 
drogenase. 



^ The conversion of glucose into ethanol is called alcoholic fermentation. 
jThe net result of this anaerobic process is 

glucose + 2 Pi + 2 ADP + 2 H H 



2 ethanol + 2 C0 2 + 2 ATP + 2 H 2 0 

It is important to note that NAD + and NADH do not appear in this 
equation, even though they are crucial for the overall reaction. NAD + 
generated in the reduction of acetaldehyde to ethanol is consumed in the 
Sxidation of glyceraldehyde 3-phosphate. Thus, there is no net oxidation- 
'ieduction in the conversion of glucose into ethanol 

-% Lactate is normally formed from pyruvate in a variety of microorgan- 
isms. The reaction also occurs in the cells of higher organisms when the 
sffiriount of oxygen is limiting, as in muscle during intense activity. The 
Induction of pyruvate by NADH to form lactate is catalyzed by lactate 
^ydrogenase. 
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I'he overall reaction in the conversion of glucose into lactate is 
Glucose + 2 Pi + 2 ADP > 2 lactate + 2 ATP + 2 H 2 0 
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Nicotinamide- 
binding half 




Adenine- 
binding half 

Figure 19-12 

Schematic diagram of the NAD + - 
binding region in dehydrogenases. 
The nicotinamide-binding half 
(shaded green) is structurally similar 
to the adenine-binding half (shaded 
yellow) of the site. [After M.G. Ross- 
mann, A. Liljas, C.-I. Branden, and 
LJ. Banaszak. In The Enzymes, vol. 10, 
3rd ed. (Academic Press, 1975), p. 68.] 



As in alcohohc fermentation, there is no net oxidation-reduction The 
NADH formed in the oxidation of glyceraldehyde 3-phosphate is con 
sumed in the reduction of pyruvate. The regeneration ofNAD + in the reduc 
Hon of pyruvate to lactate or ethanol sustains the continued operation of glycolysis 
under anaerobic conditions. If NAD+ were not regenerated, glycolysis could 
not proceed beyond glyceraldehyde 3-phosphate ( which means that no 
ATP would be generated. In effect, the formation of lactate by aerobic 
organisms buys time, as will be discussed in Chapter 22. 

3. Only a small fraction of the energy of glucose is released in its anaer- 
obic conversion into lactate (or ethanol). Much more energy can be ex 
tracted aerobically by means of the citric acid cycle and the electron trans- 
port chain. The entry point to this oxidative pathway is acetyl coenzyme A 
(acetyl CoA), which is formed inside mitochondria by the oxidative decar- 
boxylation of pyruvate. 

Pyruvate + NAD + + CoA » acetyl CoA + C0 2 + NADH 

This reaction, which is catalyzed by the pyruvate dehydrogenase complex 
will be discussed in detail in the next chapter. The NAD + required for 
this reaction and for the oxidation of glyceraldehyde 3-phosphate is re- 
generated when NADH ultimately transfers its electrons to 0 2 through 
the electron transport chain in mitochondria. 



THE BINDING SITE FOR NAD 4 " IS VERY SIMILAR 
IN MANY DEHYDROGENASES 

Lactate dehydrogenase from skeletal muscle, a 140-kd tetramer, and alco- 
hol dehydrogenase, an 84-kd dimer, have quite different three-dimen- 
sional structures. However, their NAD + -binding domains are strikingly 
simdar (Figure 19-12). This nucleotide-binding region is made up of four 
or helices and a sheet of six parallel /3 strands. The conformations of 
NAD bound to lactate dehydrogenase and to alcohol dehydrogenase 
also are nearly the same. The adenine moiety of NAD+ is bound in a 
hydrophobic crevice. In contrast, the nicotinamide unit is bound so that 
the reactive side of the ring is in a polar environment, whereas the other 
side makes contact with hydrophobic residues of the enzyme. The bound 
NAD has an extended conformation (Figure 19-13). The three-dimen- 



Figure 19-13 

Model of NAD+ The conformation 
shown here is the one found in the 
complex of NAD + and lactate dehy- 
drogenase. 




Pyrophosphate 
unit 



sional structures of two other NAD + -requiring enzymes — glycer aldehyde 
3-phosphate dehydrogenase and malate dehydrogenase (an enzyme of 
the citric acid cycle, p. 512) — are also known at high resolution. Their 
NAD """-binding domains closely resemble those of lactate dehydrogenase 
and alcohol dehydrogenase. The NAD + -binding region common to these four 
enzymes is a fundamental structural motif of NAD + -linked dehydrogenases. 



INDUCED FIT IN HEXOKINASE: 

GLUCOSE CLOSES THE ACTIVE-SITE CLEFT 

X-ray crystallographic studies of yeast hexokinase have revealed that the 
binding of glucose leads to a large conformational change in the enzyme. 
Hexokinase consists of two lobes, which come together when glucose is 
bound (Figure 19-14). Glucose induces a 12-degree rotation of one lobe 
r with respect to the other, resulting in movements of the polypeptide back- 
- bone of as much as 8 A. The cleft between the lobes closes, and the bound 
: glucose becomes surrounded by protein, except for its 5-hydroxymethyl 
y, group. 

■ \i The closing of the cleft in hexokinase is a striking example of the role 
;^ of induced fit in enzyme action, as originally proposed by Koshland 
i|; (p. 191). The glucose-induced structural changes are likely to be signifi- 
,|; cant in two ways. First, the environment around the glucose becomes 
v ,much more nonpolar, which encourages the donation of the terminal 
phosphoryl group of ATP. Second, the embracing of glucose by hexoki- 
f nase enables the enzyme to discriminate against H2O as a substrate. If 
|hexokinase were rigid, a water molecule occupying the binding site for 
the -CH 2 OH of glucose would attack the y-phosphoryl group of ATP. In 
Mother words, a rigid kinase would necessarily be an ATPase as well as a 
t kinase. The undesirable ATPase activity is prevented by making hexoki- 
nase enzymatically active only when glucose closes the cleft. It is interest- 
s' ing to note that pyruvate kinase, phosphoglycerate kinase, and phospho- 
fructokinase also contain clefts between lobes that close when substrate is 
: bound. Substrate -induced cleft closing is likely to be a general feature of kinases. 

ALDOLASE FORMS A SCHIFF BASE 
JyVITH DIHYDROXYACETONE PHOSPHATE 

Let us turn now to the catalytic mechanism of aldolase. For convenience, 
^e will examine the condensation of dihydroxyacetone phosphate and 
%lyceraldehyde 3-phosphate to form fructose 1,6-bisphosphate, the re- 
verse of the glycolytic reaction. First, dihydroxyacetone phosphate forms 
a protonated Schiff base with a specific lysine residue in the active site of 
animal aldolases. 

CH 2 OP0 3 2 - CH 2 OP0 3 2 " 
E — NH 2 + 0=C + H + ;==± E — N=C + H 2 0 

CH 2 OH H CH 2 OH 

Protonated 
Schiff base 

n this reaction, a nucleophile (the amino group) attacks the carbonyl 
roup to form a tetrahedral intermediate, which then dehydrates. The 
Tesulting protonated Schiff base plays a critical role in catalysis. It promotes the 
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Figure 19-14 

The conformation of hexokinase 
changes markedly on binding glucose 
(shown in red) . The two lobes of the 
enzyme come together and surround 
the substrate. [Courtesy of 
Dr. Thomas Steitz.] 
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formation of the enolate anion of dihydroxyacetone phosphate by serving as 
electron sink (a potent electron acceptor). 



H+ CH 2 OPCV- h+ CH 2 OP(V- 

E— N=C E-N=C +H + 

H-C-OH H-C-OH 

I e 

H 

Enolate anion 

The enolate anion then adds to the aldehyde group of glyceraldehyde 
3-phosphate to form a protonated Schiff base of fructose 1,6-bisphos- 
phate. 
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The Schiff base is deprotonated and hydrolyzed to yield fructose 1,6- 
bisphosphate and the regenerated enzyme. 

H + CH 2 0P0 3 *- H+ CH 2 OP0 3 2- 
E-N-C -^U E -N=C CH 2 0P0 3 2 - + E— NH 2 

R R Hl O 0=C 

p 

H OH OH 

I I I 
r = — C C — C— CH 2 OP0 3 2 " 

OH H H 

The pathway for the cleavage of fructose 1,6-bisphosphate is simply the 
reverse of the one for its formation. 



KINETIC PERFECTION IN CATALYSIS: 
TRIOSE PHOSPHATE ISOMERASE IN ACTION 



Figure 19-15 

Catalytic mechanism of triosephos- 
phate isomerase. 
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Much is known about the catalytic mechanism of triosephosphate isomer- 
ase (TIM), an especially interesting enzyme. TIM catalyzes the transfer of 
a hydrogen atom from C-l to C-2 in converting dihydroxyacetone phos- 
phate into glyceraldehyde 3-phosphate, an intramolecular oxidation- 
reduction. This isomerization of a ketose into an aldose proceeds through 
an enediol intermediate (Figure 19-15). When the reaction is carried out in 
D 2 0, deuterium becomes incorporated into C-2. This finding rules out a 
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direct transfer of a hydride ion (:H") from C-l to C-2. Rather, the first 
step is the removal of a proton from C-l by a basic group of the protein to 
form an enediol. Isotope labeling studies of the reverse reaction provided 
another key clue. Some of the tritium attached to C-2 of specif ically la- 
beled glyceraldehyde 3-phosphate emerged in C-l of dihydroxyacetone 
phosphate, whereas the rest exchanged with protons of water. This find- 
ing revealed that abstraction and donation of a proton are mediated by the same 
catalytic base. 

X-ray crystallographic studies then showed that glutamate 165 plays 
this role (Figure 19-16). However, a carboxylate group by itself is not basic 
enough to pull a proton away from a carbon atom adjacent to a carbonyl 
group. Histidine 95 assists catalysis by donating a proton to the C-2 car- 
bonyl group. The phosphate group of the substrate is held in place by a 
salt bridge with lysine 12 and hydrogen bonds with two main-chain NH 
groups. The positively charged end of an a helix dipole is also trained on 
the dianionic phosphate moiety. 
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Glu 165 




His 95 



PGH 



Lys 12 



Figure 19-16 

Active site of triosephosphate isomer- 
ase. Phosphoglycolhydroxamate 
(PGH), a substrate analog, is tightly 
bound. Glu 165 and His 95 play key 
catalytic roles. Lys 12 forms a salt 
bridge with the phosphate group of 
the substrate analog and polarizes the 
C-2 oxygen atom. [Drawn from 
7tim.pdb. R.C Davenport, P.A. Bash, 
B.A. Seaton, M. Karplus, G.A. Petsko, 
and D. Ringe. Biochemistry 
30(1991):582L] 



' Two other features of this isomerase are noteworthy. First, TIM displays 
• great catalytic prowess. It accelerates isomerization by a factor of 10 10 
! fbmpared with that obtained with a simple base catalyst such as acetate 
|4on. Indeed, the ^at/^M ratio for isomerization of glyceraldehyde 
;|3-phosphate is 2 x 10 8 m" 1 s" 1 , which is close to the diffusion-controlled 
pmit. In other words, the rate-limiting step in catalysis is the diffusion- 
g^pntrolled encounter of substrate and enzyme. TIM is an example of a 
mneticaUy perfect enzyme (p. 195). Second, TIM suppresses an undesired 
pde reaction, the decomposition of the enediol intermediate into 
pethylglyoxal and ? { (Figure 19-17). In solution, this unfavorable reac- 
|p£n is 100 times faster than isomerization. Hence, TIM must prevent the 
||nediol from leaving the enzyme. This labile intermediate is sequestered 
I|.the active site by the movement of a loop of 10 residues. This lid shuts 
pie active site when the enediol is present, and reopens when isomeriza- 
p£ q is com P leted - We see here a striking example of the channeling of catalysis by- 
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Figure 19-17 

In solution, the enediol intermediate 
in the isomerization of dihydroxyace- 
tone phosphate decomposes into 
methylglyoxal and Pi. The isomerase 
blocks this undesired side reaction. 



THIOESTER IS FORMED IN THE OXIDATION 
~ GLYCERALDEHYDE 3-PHOSPHATE 



j^eraldehyde 3-phosphate dehydrogenase catalyzes the oxidative phos- 
Qorylation of its aldehyde substrate. 



lyceraldehyde 3-phosphate + ? { + NAD+ 



1,3-BPG + NADH + H 
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The conversion of an aldehyde into an acyl phosphate entails both „„v 

»on^ phos p horyla « Qxidation reqU iL^e p !^:r 0 /:^t: 

nviide f 15 " hyd ,T gen nUClCUS 3nd °"> electrons/Removal 
hydnde ton from an aldehyde is energetically costly because of the di 
polar character of the carbonyl group: the carbon aim of the ca/bon 
group already has a partial positive charge. Removal of the hydride Ton 
gready facilitated by making the carbon atom less positively chafed Th 
1S accomphshed by the addition of a nucleophile, represented he 2 by 



H + + 0=C-H + x- 
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R 

0=C + : H~ + H 4 
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The hydride ion readily leaves the addition compound because the 
carbon atom no longer carries a large positive charge. FuXmore 

IT ahh C free f e ^ ° f ^ ° Xidadon is P— - inAe acyfin ™ ' 
diate. Addition of orthophosphate to this acyl intermediate yields an^vl 
phosphate, which has a high group transfer potential Y Y 
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Now let us see how glyceraldehyde 3-phosphate dehydrogenase carries 
out these steps (Figure 19-18). The nucleophik X~ isthesulZd^Zofl 
cystine residue ■at the active site. The aldehyde substrate ^l ^^. 
zed form of this sulfhydryl group to form a hemithioacetal. The nZ Z 

enzyme. I he products of this reaction are the reduced coenzyme NADH 

the acyl intermediate mentioned earlier. NADH dissociates from fh/™ 
zyme and another NAD + binds to the active die^SSSS^S 
attacks the thioester to form 1,3-BPG, a high-potential pho^phoryl donor 




CONH, 



Figure 19-18 

Catalytic mechanism of 
glyceraldehyde 3-phosphate 
dehydrogenase. 




CONK 



CONH, 



A crucial aspect of the formation of 1,3-BPG from glyceraldehyde 
3-phosphate is that a thermodynamically unfavorable reaction, the forma- 
tion of an acyl phosphate from a carboxylate, is driven by a thermody- 
namically favorable reaction, the oxidation of an aldehyde 
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? ° 
R-C-H + NAD + + H 2 0 ;=± R-C-CT + NADH + 2 H+ 

These two reactions are coupled by the thioester intermediate, which preserves 
much of the free energy released in the oxidation reaction. We see here 
the use of a covalent enzyme-bound intermediate as a mechanism of energy cou- 
pling. J ° J 

ARSENATE, AIM ANALOG OF PHOSPHATE, POISONS BY 
UNCOUPLING OXIDATION AND PHOSPHORYLATION 

Arsenate (As0 4 3 -) closely resembles P £ in structure and reactivity In the 
reaction catalyzed by glyceraldehyde 3-phos P hate dehydrogenase, arse- 
nate can replace phosphate in attacking the energy-rich thioester inter- 
mediate. The product of this reaction, l-arseno-3-phosphoglycerate is 
unstable, m contrast with 1,3-bisphosphoglycerate. l-Arseno-3-phospho- 
glycerate and other acyl arsenates are rapidly and spontaneously hydro- 
lyzed. Hence, the net reaction in the presence of arsenate is 
Glyceraldehyde 3-phosphate + NAD+ + H 2 0 



3-phosphoglycerate + NADH + 2 H + 
Note that glycolysis proceeds in the presence of arsenate but that the 
ATP normally formed in the conversion of 1,3-bisphosphoglycerate into 
3-phosphoglycerate is lost. Thus, arsenate uncouples oxidation and phosphory- 
lation by forming a highly labile acyl arsenate. Arsenic is a potent poison 
because arsenate generally substitutes for phosphate in phosphoryl trans- 
fer reactions. Also, arsenite (AsO-T) forms adducto with thiols. One likely 
reason for the choice of phosphorus over arsenic in the evolution of bio- 
molecules is the greater kinetic stability of its energy-rich compounds 
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H— C— OH 

CH 2 OP0 3 2 - 
l-Arseno-3-phosphoglycerate 



2.3-BISPHOSPHOGLYCERATE, AN ALLOSTERIC EFFECTOR OF 
HEMOGLOBIN, ARISES FROM 1,3-BISPHOSPHOGLYCERATE 

We have seen that fructose 2,6-bisphosphate, a regulatory molecule, 
arises from a glycolytic intermediate. Another regulatory molecule com- 
>ng irom this pathway is 2, 3 -bisphosphoglycerate (2,3-BPG), a controller of 
oxygen transport in erythrocytes (p. 160). Red blood cells have a high 
concentration of 2,3-BPG, typically 4 m M) in contrast with most other 
ceils which have only trace amounts. The synthesis and degradation of 
<>o-crG are a short detour from the glycolytic pathway (Figure 19-19). 
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Figure 19-19 

Synthesis and deg- 
radation of 2,3- 
bisphosphoglycerate. 
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Table 19-3 

Typical concentrations of glycolytic 
intermediates in erythrocytes 



Intermediate 




Glucose 


5000 


Glucose 6-phosphate 


83 


Fructose 6-phosphate 


14 


Fructose 1,6-bisphosphate 


31 


Dihydroxyacetone phosphate 


138 


Glyceraldehyde 3-phosphate 


19 


1,3-Bisphosphoglycerate 


1 


2,3-Bisphosphoglycerate 


4000 


3-Phosphoglycerate 


118 


2-Phosphoglycerate 


30 


Phosphoenolpyruvate 


23 


Pyruvate 


51 


Lactate 


2900 


ATP 


1850 


ADP 


138 


P f 


1000 



After S. Minakami and H. Yoshikawa. 
Biochem. Biophys. Res. Comm. 18(1965):345. 



Bisphosphoglycerate mutase converts 1,3-BPG into 2,3-BPG. 2,3-BPG is a po- 
tent competitive inhibitor of its own formation. The concentration of this 
allosteric effector also depends on its hydrolysis to 3-phosphoglycerate, a 
reaction catalyzed by 2 \ 3 -bisphosphogly cerate phosphatase. Both the synthesis 
and degradation of 2,3-BPG are nearly irreversible. 

The mutase reaction has an interesting mechanism. 3 -Phosphogly cerate is 
an obligatory participant although it does not appear in the overall stoi- 
chiometry. The mutase binds 1,3-BPG and 3-phosphoglycerate simulta- 
neously. In this ternary complex, a phosphoryl group is transferred from 
position 1 of 1,3-BPG to position 2 of 3-phosphoglycerate (Figure 19-20). 



+ 



+ 



a-® 




Figure 19-20 

3-Phosphoglycerate participates in the conversion of 1,3-bisphosphoglycerate into 
2,3-bisphosphoglycerate. 



2,3-BPG is important not only in red blood cells. Phosphoglycerate* 
mutase, the enzyme catalyzing the interconversion of 3-phosphoglycerate 
and 2-phosphoglycerate, requires a catalytic amount of 2,3-BPG to be 
active. 2,3-BPG phosphorylates a histidine residue at the catalytic site to 
maintain the active form of the enzyme. Thus, 2,3-bisphosphoglycerate 
plays an essential role in basic metabolism. It seems likely that 2,3-BPG 
was present eons before it was recruited by red cells to control the oxygen 
affinity of hemoglobin. 



ENOL PHOSPHATES ARE POTENT PHOSPHORYL DONORS 

Because it is an acyl phosphate, 1,3-BPG has a high group transfer poten- 
tial. A different kind of high-energy phosphate compound is generated 
several steps later in glycolysis. Phosphoenolpyruvate, an enol phosphate, is 
formed by the dehydration of 2-phosphoglycerate. The AG 0/ of hydrolysis 
of a phosphate ester of an ordinary alcohol is -3 kcal/mol, whereas that 
of phosphoenolpyruvate is -14.8 kcal/mol. Why does phosphoenolpyru- 
vate have such a high phosphoryl group-transfer potential? The answer 
is that the enol formed upon transfer of the phosphoryl group undergoes 
a conversion into a ketone — namely, pyruvate. 



V°~ o V" V" 



C-O-P-O- _ C-OH — ^0 

CH 2 CT CH 2 CH 3 

Phosphoenolpyruvate Enolpyruvate Pyruvate 



The AG°' of the enol-ketone conversion is very large, of the order of 
-10 kcal/mol. Thus, the high phosphoryl group-transfer potential of phospho- 
enolpyruvate arises primarily from the large driving force of the subsequent enol- 
ketone conversion. 



A FAMILY OF TRANSPORTERS ENABLES GLUCOSE 
TO ENTER AND LEAVE ANIMAL CELLS 
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The thermodynamically downhill movement of glucose across the plasma 
membrane of animal cells is mediated by several glucose transporters. The 
members of this protein family, named GLUT1 to 5, consist of a single 
polypeptide chain about 500 residues long. The common structural 
theme is the presence of 12 transmembrane segments (Figure 19-21). 
The binding site for glucose alternately faces the inside and outside of the 
cell when occupied by a sugar; this eversion is accomplished by conforma- 
tional changes within the transporter and not by a rotation of the whole 
protein. 




/"The members of this family have distinctive roles: 

i \1. GLUT1 and 3, present in nearly all mammalian cells, are responsible 
rfpr basal glucose uptake. Their Km for glucose is about 1 mM, significantly 

less than the normal serum glucose level, which typically ranges from 
• 4 him to 8 mM. Hence, GLUTl and 3 continually transport glucose at an 
"essentially constant rate. 

^ ;2. GLUT5, present in the small intestine, works in tandem with the 
Na^-glucose symporter in the absorption of glucose from the gut. The 
Jsymporter pumps glucose into the intestinal epithelial cell. GLUT5 in the 
Hiasma membrane on the opposite side of the cell then releases glucose 
into the bloodstream. 

tl^3. GLUT2, present in liver and pancreatic /? cells, is distinctive in hav- 
ing a very high K M for glucose (15-20 mM). Hence, the rate of entry of 
|lucose into these tissues is proportional to the blood glucose level. The 
pancreas can thereby sense the glucose level and accordingly adjust the 
rate of insulin secretion. The high /£ M of GLUT2 also assures that glucose 
fapidly enters liver cells only in times of plenty. When the blood glucose 
||vel is low, glucose preferentially enters brain and other tissues because 
ffteir glucose transporters have a lower than that of liver. 

? 4. GLUT4, which has a K M of 5 mM, mediates the entry of glucose into 
jnuscle and fat cells. Insulin, which signals the fed state, leads to a rapid 
increase in the number of GLUT4 transporters in the plasma membrane. 
Hence, insulin promotes the uptake of glucose by muscle and fat. 

Shis family of transporters vividly illustrates how isoforms of a single pro- 
gin can profoundly shape the metabolic character of cells and contribute 
& their diversity and functional specialization. 



Figure 19-21 

Model of a mammalian glucose trans- 
porter. The hydrophobicity profile of 
the protein points to the presence of 
12 transmembrane a helices. Chemi- 
cal labeling and try p tic cleavage stud- 
ies support this postulated topography. 
[From M. Muekler, C. Caruso, S.A. 
Baldwin, M Panico, M. Blench, H.R. 
Morris, W.J. Allard, G.E. Lienhard, 
and H.F. Lodish. Science 
229(1985):94L] 



SUMMARY 



Glycolysis is the set of reactions that converts glucose into pyruvate. In 
aerobic organisms, glycolysis is the prelude to the citric acid cycle and the 
electron transport chain, where most of the free energy in glucose is 
harvested. The 10 reactions of glycolysis occur in the cytosol. In the first 
stage, glucose is converted into fructose 1,6-bisphosphate by a phosphory- 
lation, an isomerization, and a second phosphorylation reaction. Two 
molecules of ATP are consumed per molecule of glucose in these reac- 
tions, which are the prelude to the net synthesis of ATP. In the second 
stage, fructose 1,6-bisphosphate is cleaved by aldolase into dihydroxyace- 
tone phosphate and glyceraldehyde 3-phosphate, which are readily inter- 
convertible. Glyceraldehyde 3-phosphate is then oxidized and phosphor- 
ylated to form 1,3-bisphosphoglycerate, an acyl phosphate with a high 
phosphoryl transfer potential. 3-Phosphoglycerate is then formed as an 
ATP is generated. In the last stage of glycolysis, phosphoenolpyruvate, a 
second intermediate with a high phosphoryl transfer potential, is formed 
by a phosphoryl shift and a dehydration. Another ATP is generated as 
phosphoenolpyruvate is converted into pyruvate. There is a net gain of 
two molecules of ATP in the formation of two molecules of pyruvate from 
one molecule of glucose. 

The electron acceptor in the oxidation of glyceraldehyde 3-phosphate 
is NAD + , which must be regenerated for glycolysis to continue. In aerobic 
organisms, the NADH formed in glycolysis transfers its electrons to O2 
through the electron transport chain, which thereby regenerates NAD + . 
Under anaerobic conditions, NAD + is regenerated by the reduction of 
pyruvate to lactate. In some microorganisms, NAD + is normally regener- 
ated by the synthesis of lactate or ethanol from pyruvate. These two pro- 
cesses are examples of fermentations. 

The glycolytic pathway has a dual role: it degrades glucose to generate 
ATP, and it provides building blocks for the synthesis of cellular compo- 
nents. The rate of conversion of glucose into pyruvate is regulated to 
meet these two major cellular needs. Under physiologic conditions, the 
reactions of glycolysis are readily reversible except for the ones catalyzed 
by hexo kinase, phosphofructokinase, and pyruvate kinase. Phosphofruc- 
tokinase, the most important control element in glycolysis, is inhibited by 
high levels of ATP and citrate, and it is activated by AMP and fructose 
2,6-bisphosphate. In liver, this bisphosphate signals that glucose is abun- 
dant. Hence, phosphofructokinase is active when either energy or build- 
ing blocks are needed. Hexokinase is inhibited by glucose 6-phosphate, 
which accumulates when phosphofructokinase is inactive. Pyruvate ki- 
nase, the other control site, is alios terically inhibited by ATP and alanine, 
and it is activated by fructose 1,6-bisphosphate. Consequently, pyruvate 
kinase is maximally active when the energy charge is low and glycolytic 
intermediates accumulate. Pyruvate kinase, like the bifunctional enzyme 
controlling the level of fructose 2,6-bisphosphate, is regulated by phos- 
phorylation. A low level of glucose in the blood promotes the phosphory- 
lation of liver pyruvate kinase, which diminishes its activity and thereby 
decreases glucose consumption in the liver. 

Induced fit enhances the catalytic efficiency of hexokinase and triose- 
phosphate isomerase by preventing undesired side reactions. The isomer- 
ase is an example of a catalytically perfect enzyme limited only by the 
diffusion-controlled encounter of enzyme and substrate. A thioester in- 
termediate conserves some of the free energy of oxidation of glyceralde- 
hyde 3-phosphate; attack by P[ yields an energy-rich acyl phosphate. Arse- 
nate, an analog of phosphate, uncouples oxidation and phosphorylation. 
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mediated by a family of glucose transporters named GLUT1 to 5. The 

binding site for glucose in this plasma membrane protein, when occu- GLYCOLYSIS 
pied, alternately faces the inside and outside of the cell. The different 
values and differential regulation of this family of transporters shape the 
metabolic character of cells in different organs. 
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PROBLEMS 



1. Kitchen chemistry. Sucrose is commonly used to preserve 
fruits. Why is glucose not well suited for preserving 
foods? 

2. Tracing carbon atoms. Glucose labeled with 14 C at C-l is 
incubated with the glycolytic enzymes and necessary 
cofactors. 

(a) What is the distribution of 14 C in the pyruvate that 
is formed? (Assume that the interconversion of 
glyceraldehyde 3-phosphate and dihydroxyace- 
tone phosphate is very rapid compared with the 
subsequent step.) 

(b) If the specific activity of the glucose substrate is 
10 mCi/mM, what is the specific activity of the py- 
ruvate that is formed? 

3. Lactic fermentation. Write a balanced equation for the 
conversion of glucose into lactate. 

(a) Calculate the standard free -energy change of this 
reaction using the data given in Table 19-2 
(p. 491) and the fact that AG°' is -6 kcal for the 
reaction 

Pyruvate + NADH + H+ ;=± lactate + NAD + 

(b) What is the free-energy change (AG', not AG°') of 
this reaction when the concentrations of reactants 
are: glucose, 5 him; lactate, 0.05 him; ATP, 2 itim; 
ADP, 0.2 itim; and Pi, 1 mM? 

4. High potential. What is the equilibrium ratio of phos- 
phoenolpyruvate to pyruvate under standard condi- 
tions when [ATP] /[ADP] = 10? 

5. Hexose-triose equilibrium. What are the equilibrium con- 
centrations of fructose 1,6-bisphosphate, dihydroxyace- 
tone phosphate, and glyceraldehyde 3-phosphate when 
1 mM fructose 1,6-bisphosphate is incubated with aldol- 
ase under standard conditions? 

6. Double labeling. 3-Phosphoglycerate labeled uniformly 
with 14 G is incubated with 1,3-BPG labeled with 32 P at 
C-l. What is the radioisotope distribution of the 2,3- 
BPG that is formed on addition of BPG mutase? 

7. An informative analog. Xylose has the same structure as 
glucose except that it has a hydrogen atom at C-6 in 



place of a hydroxymethyl group. The rate of ATP hy- 
drolysis by hexokinase is markedly enhanced by the 
addition of xylose. Why? 

8. The far reach of glycolysis. Oxygen transport can be af- 
fected in genetic disorders of glycolysis in red cells. 

(a) How are glycolysis and oxygen transport linked? 

(b) How is oxygen affinity altered by a deficiency of 
hexokinase? 

(c) How is oxygen affinity altered by a deficiency of 
pyruvate kinase? 

9. Distinctive sugars. The intravenous infusion of fructose 
into healthy volunteers leads to a two- to fivefold in- 
crease in the level of lactate in the blood, a far greater 
increase than that observed following the infusion of 
the same amount of glucose. 

(a) Why is glycolysis more rapid following the infusion 
of fructose? 

(b) Fructose has been used in place of glucose for in- 
travenous feeding. Why is this use of fructose un- 
wise? 

10. Catalytic metal. Aldolases in prokaryotes contain a 
tightly bound divalent metal ion that is essential for 
catalysis. Propose a catalytic function for this metal ion. 

11. Contrasting inactivators. Prokaryotic aldolases are inacti- 
vated by ethylenediaminetetraacetate (EDTA), a chela- 
tor of divalent metal ions, whereas animal aldolases dis- 
cussed in this chapter are inactivated by sodium 
borohydride. Account for this difference. 

12. Metabolic mutants. Predict the effect of each of the fol- 
lowing on the pace of glycolysis in liver cells: 

(a) Loss of the allosteric site for ATP in phosphofructo- 
kinase. 

(b) Loss of the binding site for citrate in phosphofruc- 
tokinase. 

(c) Loss of the phosphatase domain of the bifunc- 
tional enzyme that controls the level of fructose 
2,6-bisphosphate. 

(d) Loss of the binding site for fructose 1,6-bisphos- 
phate in pyruvate kinase. 



H 2 N — C=NH 2 + 
I 

NH 
I 

CH 2 

CH 2 

CH 2 

H — C — NH 3 + 
I 

coo- 

Arginine 



H 2 0 C 
H.N— C 



-NH 2 + 



Urea 



NH 3 + 

CH 2 
I 2 
CH 2 

CH 2 

I 

H— C— NH + 

I 

coo- 

Ornithine 



The other reactions of the urea cycle lead to the synthesis, of arginine 
from ornithine. First, a carbamoyl group is transferred to ornithine to 
form citrulline, in a reaction catalyzed by ornithine transcarbamoylase. The 
carbamoyl donor in this reaction is carbamoyl phosphate, which has a high 
transfer potential because of its anhydride bond. 
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^Argininosuccinate synthetase then catalyzes the condensation of citrulline 
;ai*a aspartate. This synthesis of argininosuccinate is driven by the cleavage 
of ATP into AMP and pyrophosphate and by the subsequent hydrolysis of 
.pyrophosphate. 
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ipally, argininosuccinase cleaves argininosuccinate into arginine and fuma- 
J?. These two reactions, which transfer the amino group of aspartate to 
fepn arginine, preserve the carbon skeleton of aspartate in the form of 
tLmarate. 
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Figure 28-7 

Biosynthetic families of amino acids 
in bacteria and plants. Major meta- 
bolic precursors are shaded blue. 
Amino acids that give rise to other 
amino acids are shaded red. Essential 
amino acids are shown in boldface. 
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The nonessential amino acids are synthesized by quite simple reac- 
tions whereas the pathways for the formation of the essential amino acids 
are quite complex. For example, the nonessential amino acids alanine and 
aspartate are synthesized in a single step from pyruvate and oxaloacetate, 
respectively. Each acquires its amino group from glutamate in a transami- 
nation reaction in which pyridoxal phosphate is the cofactor (p. 631): 

± alanine + a-ketoglutarate 
± aspartate + a-ketoglutarate 



Pyruvate + glutamate 
Oxaloacetate + glutamate 



Asparagine is then synthesized by the amidation of aspartate: 

Aspartate + NH 4 + + ATP > asparagine + AMP + PPi + H + 

In mammals, the nitrogen donor in the synthesis of asparagine is gluta- 
mine rather than NH 4 + , as in bacteria. Ammonia generated at the active 
site of the enzyme is directly transferred to bound aspartate. Recall that 
high levels of NH 4 + are toxic to humans (p. 637) . 

Another one-step synthesis of a nonessential amino acid in mammals is 
the hydroxylation of phenylalanine (an essential amino acid) to tyrosine. 

Phenylalanine + 0 2 + NADPH + H + > tyrosine + NADP + + H 2 0 

This reaction is catalyzed by phenylalanine hydroxylase, a monooxygenase 
discussed previously (p. 647). It is noteworthy that tyrosine is an essenua. 
amino acid in people lacking this enzyme. 



GLUTAMATE IS THE PRECURSOR OF GLUTAMINE, 
PROLINE, AND ARGININE 

The synthesis of glutamate by the reductive amination of a-ketoglutarate 
has already been discussed, as has the conversion of glutamate into gtttf 
mine (p 717). Glutamate is the precursor of two other nonessential ami' 
acids,. proline and arginine. First, the ^carboxyl group of glutamate reac 
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